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ABSTRACT 


The  CMARC  panel  code  is  evaluated  to  verify  its  accuracy  and  suitability  for  the 
development  of  an  aerodynamic  model  of  the  Naval  Postgraduate  School  (NPS)  FROG 
Unmanned  Air  Vehicle  (UAV).  CMARC  is  a  DOS  personal  computer  based  version  of 
the  NASA  Panel  Method  Ames  Research  Center  (PMARC)  panel  code.  The  core 
processing  algorithms  in  CMARC  are  equivalent  to  PMARC.  CMARC  enhancements 
include  improved  memory  management  and  command  line  functionality.  Both  panel  codes 
solve  for  inviscid,  incompressible  flow  over  complex  three-dimensional  bodies  using 
potential  flow  theory.  Emphasis  is  first  placed  on  verifying  CMARC  against  the  PMARC 
and  NPS  Unsteady  Potential  Flow  (UPOT)  panel  codes.  CMARC  boundary  layer 
calculations  are  then  compared  to  experimental  data  for  an  inclined  prolate  spheroid. 
Finally,  a  complex  three-dimensional  panel  model  is  developed  for  aerodynamic  modeling 
of  the  FROG  UAV.  CMARC  off-body  flow  field  calculations  are  used  to  generate  static- 
source  and  angle-of-attack  vane  position  corrections.  Position  corrections  are  provided  in 
look-up  table  and  curve  fit  formats.  Basic  longitudinal  and  lateral-directional  stability 
derivatives  are  also  developed  with  CMARC  data.  CMARC  derived  stability  derivatives 
are  sufficiently  accurate  for  incorporation  into  an  initial  aerodynamic  model.  Adjustments 
through  analysis  of  flight  test  data  may  be  required.  Future  CMARC  studies  should 
concentrate  on  the  development  of  the  damping  and  control  power  derivatives. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Computational  fluid  dynamics  (CFD)  is  increasingly  used  as  a  design  and  analysis 
tool.  As  the  price  of  computer  hardware  drops  and  computational  power  increases,  CFD 
becomes  more  attractive  to  a  larger  audience.  CFD  tools  range  from  the  high  end  three- 
dimensional  (3D)  Navier-Stokes  solvers  for  compressible,  viscous  fluids  to  potential  flow 
solvers  for  incompressible,  inviscid  flows.  This  paper  discusses  the  development  of  a  DOS 
personal  computer  hosted  panel  code  model  for  the  Naval  Postgraduate  School  (NFS) 
Fiber  Optic  Guided  (FROG)  Unmanned  Air  Vehicle  (UAV)  program. 

The  Personal  Simulation  Works  software  suite,  consisting  of  LOFTSMAN, 
CMARC  and  POSTMARC,  is  used  for  all  aspects  of  the  study.  The  software  provides  for 
panel  model  development,  input  file  processing  and  the  visualization  of  results.  Emphasis 
is  placed  on  verifying  both  the  accuracy  and  suitability  of  the  CFD  programs  for 
aerod5mamic  modeling. 

Until  recently,  personal  computers  (PC)  did  not  have  the  computational  power  or 
memory  to  be  practical  for  panel  code  CFD  programs.  Things  have  changed  with  the 
introduction  of  the  Pentium  class  PC  and  low  cost  RAM.  AeroLogic  capitalized  on  the 
power  of  the  Pentium  class  PC  and  developed  Personal  Simulation  Works  (PSW).  PSW  is 
centered  around  the  3D  low  order,  inviscid  potential  flow  solver  named  CMARC, 
CMARC  is  a  re-hosted  version  of  NASA’s  Panel  Method  Ames  Research  Code 
(PMARC).  PMARC  was  re-written  in  the  C  language  and  compiled  for  IBM  compatible 
PCs.  CMARC  runs  under  the  DOS  operating  system.  CMARC  will  also  run  in  a  DOS 
window  under  the  WINDOWS  3.x,  95  or  NT  operating  systems.  CMARC  has  enhanced 
capabilities  that  include;  improved  memory  management,  an  expanded  set  of  command 
line  switches  and  provisions  for  expanded  boundary  layer  post-processing  capabilities. 
However,  the  core  processing  algorithms  remain  the  same  as  implemented  in  PMARC. 

LOFTSMAN,  the  PSW  pre-processing  program,  is  used  to  mesh  complex  3D 
bodies  and  create  input  file  patches.  The  program  runs  under  the  DOS  operating  system 
and  allows  the  user  to  loft  conics  based  3D  surfaces.  The  program  automatically  creates 
CMARC,  PMARC  or  VSAERO  input  patches  based  on  desired  panel  densities  and 
distribution. 
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POSTMARC  is  used  for  flow  visualization  and  integration  of  resultant  forces.  It 
runs  under  the  Windows  3.x,  95  and  NT  operating  systems.  POSTMARC  reads  CMARC 
or  PMARC  output  files  and  provides  for  the  visualization  of  model  geometry,  wake 
stepping,  on  and  off-body  streamlines  and  surface  phenomena. 

B.  REQUIREMENTS 

The  Naval  Postgraduate  School  Aeronautics  Department  is  integrating  UAV 
hardware  and  software  to  demonstrate  autonomous  flight,  trajectory  tracking  and 
automatic  landing.  A  core  requirement  for  flight  control  development  is  a  valid 
aerodynamic  truth  model  for  the  UAV  airfi-ame.  The  introduction  of  each  new  airframe 
requires  the  development  of  a  new  aerodynamic  truth  model.  Most  recently, 
Papageorgiou  [Ref  1]  developed  and  tested  an  aerodynamic  model  for  the  NPS  FROG 
UAV  based  on  classical  methods.  His  model  produced  a  close  match  to  flight  test  results 
in  the  longitudinal  axis.  However,  the  lateral-directional  axis  required  modifications  based 
on  measured  flight  test  data  to  produce  acceptable  results.  With  the  avmlability  of  low 
cost  panel  code  CFD  capabilities,  it  is  suggested  that  a  panel  code  model  of  the  FROG 
UAV  is  an  alternative  for  estimating  many  of  the  stability  derivatives  required  for  an 
aerodynamic  truth  model. 

Accurate  pitot-static  and  angle-of-attack  sensors  are  required  for  highly 
augmented  flight  control  systems.  CMARC  is  well  suited  for  solving  on-body  static 
pressure  distributions  and  off-body  flow  velocities  over  the  predominately  attached  flow 
fields  of  fuselage  fore  bodies.  This  proves  particularly  useful  for  generating  pitot  static 
and  angle  of  attack  correction  curves  and  look-up  tables. 

C.  STATEMENT  OF  OBJECTIVES 

The  Naval  Postgraduate  School  Aeronautics  Department  has  both  active  CFD 
research  and  avionics  development  programs.  The  primary  purpose  of  this  investigation  is 
to  verify  the  accuracy  and  suitability  of  the  PSW  software  suite  while  developing  a  panel 
code  model  for  the  NPS  FROG  UAV  program.  Specific  objectives  are  as  follows: 

•  Demonstrate  panel  code  modeling,  processing  and  visualization  on  a  Pentium 
PC  using  the  PSW  software  package. 
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•  Verify  CMARC  results  against  PM  ARC. 

•  Investigate  the  CMARC  integral  boundary  layer  calculations  through 
comparison  to  validated  2D  CFD  codes  and  3D  experimental  data. 

•  Develop  and  analyze  a  panel  code  model  for  the  NPS  FROG  UAV  using  PSW 
to  estimate  basic  stability  derivatives  and  produce  angle-of-attack  vane  and 
pitot-static  correction  curves. 

•  Compare  relative  speed  of  CMARC  hosted  on  1 5  0  MHz  Pentium  personal 
computer  to  PMARC  hosted  on  a  Silicon  Graphics  Indigo^  workstation. 
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n.  OVERVIEW  OF  PERSONAL  SIMULATION  WORKS 

A.  GENERAL 

Personal  Simulation  Works  is  a  PC  based  software  suite  that  provides  for  the  three 
primary  CFD  requirements;  3D  modeling  of  an  aircraft  (LOFTSMAN),  panel  code  flow 
solver  (CMARC),  and  post-processing  of  the  computed  flow  field  (POSTMARC).  The 
software  package  contains  three  applications  hosted  on  the  IBM  compatible  personal 
computer.  Each  software  program  is  discussed  separately. 

B.  LOFTSMAN 

LOFTSMAN  is  a  3D  modeling  tool  that  generates  surface  panel  distributions  for 
CMARC  or  PMARC  input  files.  The  program  is  based  on  conics,  which  allows  rapid 
lofting  of  streamlined  bodies  such  as  aircraft  fuselages  and  engine  nacelles.  In  addition, 
wing  and  control  surfaces  can  be  designed  with  the  extensive  library  of  airfoil  templates  or 
with  user  specified  coordinates.  The  software  is  well  documented,  including  a  tutorial,  in 
the  Personal  Simulation  Works  User  Guide  [Ref  2].  LOFTSMAN  is  primarily  designed 
for  creating  new  objects,  but  an  existing  airframe  can  be  matched  quite  closely  with  just  a 
detailed  three-view  drawing  that  includes  frame  cross  sections. 

1.  Streamlined  Bodies 

LOFTSMAN  functionality  is  divided  into  Body  Objects  and  Wing  Objects.  In 
general,  they  remain  separate  unless  the  intersection  between  a  wing  and  body  is  required. 

Body  Objects  are  created  using  a  family  of  curves  called  second-degree  conics. 
Circles,  ellipses,  parabolas  and  hyperbolas  are  among  this  group.  An  entire  fuselage  is 
described  by  specifying  just  four  lines.  These  are  the  top  waterline  (TW),  bottom 
waterline  (BW),  the  maximum  breadth  line  (MB)  and  the  waterline  of  the  maximum 
breadth  line  (WW).  For  each  line,  the  beginning,  ending  and  a  few  points  along  the  line 
are  specified.  Control  points  are  also  specified  with  a  curvature  factor  that  allows 
LOFTSMAN  to  generate  a  smooth  conic  between  the  points.  The  power  of  conic  lofting 
will  become  evident  when  discussing  the  modeling  of  the  complex  FROG  UAV  fuselage  in 
Chapter  V. 
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2. 


Wings  and  Control  Surfaces 


Wings  and  control  surfaces  are  easily  specified  in  LOFTSMAN  using  a  short  input 
file  created  with  any  text  editor.  The  file  specifies  root,  intermediate  and  tip  rib  section, 
location,  axis,  chord  and  incidence.  LOFTSMAN  then  fairs  a  smooth  surface  through  the 
rib  sections.  Washout  is  specified  by  varying  the  incidence  of  the  root  and  tip  ribs. 
Sweep-back  is  controlled  by  staggering  the  tip  rib  location  with  respect  to  the  root  rib. 
Once  the  general  wing  surface  is  specified,  control  surfaces  such  as  ailerons,  flaps  and 
elevators  can  be  deflected  and  meshed. 

3.  Patches 

LOFTSMAN  automatically  meshes  3D  surfaces  and  creates  patches  for 
CMARC/PMARC  input  files.  The  distinction  between  a  mesh  and  a  patch  is  important.  A 
mesh  is  a  set  of  quadrilateral  and  triangular  panels  that  represent  the  surface  of  a  wing  or 
body.  When  the  set  of  panels  is  organized  and  formatted  to  create  a  sub-component 
portion  of  a  CMARC  or  PMARC  input  file,  it  is  called  a  patch. 

A  body  or  wing  surface  is  first  meshed  at  a  density  specified  by  the  user.  Cosine 
and  half-cosine  spacing  are  among  the  compression  options.  After  meshing  the  object, 
one  saves  it  to  a  text  file  as  a  formatted  patch.  One  then  opens  the  patch  file  with  any  text 
editor  and  copies/pastes  the  patch  text  into  the  appropriate  location  in  the  CMARC  input 
file. 

Each  control  surface  deflection  requires  a  separate  mesh  and  formatted  patch.  For 
instance,  to  evaluate  roll  performance  one  needs  to  separately  mesh  an  upward  aileron 
deflection  on  the  right  wing  and  a  downward  deflection  on  the  left  wing.  If  multiple 
deflections  of  a  single  control  surface  are  required,  each  deflection  must  be  meshed 
separately. 

C.  CMARC 

CMARC  is  the  C  version  of  the  Panel  Method  NASA  Ames  Research  Center 
(PMARC)  low  order,  3D  panel  code.  Inviscid,  irrotational,  incompressible,  potential  flow 
is  assumed.  Low  order  means  that  source  and  doublet  strength  distribution  is  constant 
across  a  panel.  There  is  no  attempt  to  match  the  source  or  doublet  strength  of  an  adjacent 
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panel  at  a  common  edge.  Advanced  features  include  internal  flow  modeling  and  time 
stepping  wake  models. 

PMARC  version  12.19  was  released  as  FORTRAN  77  source  code  in  1992. 
CMARC  was  rewritten  in  the  C  language  and  compiled  for  hosting  on  IBM  compatible 
personal  computers  by  AeroLogic,  Inc.  The  program  runs  under  the  DOS  operating 
system.  It  will  also  run  in  a  DOS  window  from  Windows  3.1,  95  or  NT.  Enhanced 
features  include  command  line  options  and  flexible  memory  management.  Command  line 
options  simplify  batch  processing  by  adding  an  extensive  set  of  switches  that  can  be  set 
external  to  the  CMARC  input  file.  Flexible  memory  management  provides  for  the 
automatic  sizing  of  arrays  without  having  to  recompile  the  source  code. 

D.  POSTMARC 

POSTMARC  is  a  Windows  post-processing  program  for  the  visualization  of 
CMARC  and  PMARC  output  files.  Capabilities  include  body  geometry,  wake  stepping, 
surface  pressure  and  streamline  visualization.  POSTMARC  also  provides  the  capability  to 
integrate  pressure  and  skin  fnction  forces  over  the  model  geometry.  This  proves 
particularly  useful  when  one  desires  to  recalculate  loads  around  a  different  center  of 
gravity. 

An  interesting  feature  for  design  work  integrates  panel  surface  area  to  obtain  total 
wetted  area.  After  lofting  a  new  geometry  in  LOFTSMAN,  a  quick  check  of  geometry  is 
made  by  running  CMARC  with  the  -g  command  line  toggle.  The  total  wetted  area  is  then 
checked  in  POSTMARC.  This  function  is  particularly  useful  when  working  to  reduce  skin 
friction  drag. 

Versions  1.17.3  and  later  of  POSTMARC  include  the  capability  to  integrate  skin 
fiiction  drag  coefficient  over  the  model  geometry.  It  is  important  to  note  that  a  key  piece 
of  the  drag  equation  is  missing  from  a  POSTMARC  solution.  CMARC  provides  induced 
drag  from  the  surface  pressure  distribution  and  skin  fnction  drag  from  the  2D  boundary 
layer  code.  Skin  fnction  is  only  calculated  up  to  the  point  of  boundary  layer  separation. 
Pressure  drag  due  to  separation,  a  major  portion  of  the  drag  equation,  is  missing  from  a 
CMARC/POSTMARC  solution. 

In  fact,  if  one  isn’t  careful,  POSTMARC  drag  calculations  can  be  misleading. 
Take  for  instance  two  similar  model  configurations  with  only  minor  geometry  differences 
that  do  not  affect  wetted  area.  It  is  possible  for  the  model  with  more  flow  separation  to 
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have  less  skin  friction  drag  because  there  is  no  CMARC  output  for  skin  friction  coefficient 
after  the  boundary  layer  code  predicts  separation.  During  iterative  design  work,  this  could 
lead  to  the  incorrect  conclusion  that  the  design  team  is  reducing  overall  drag.  Perhaps  a 
better  function  for  LOFTSMAN  than  integrated  skin  fnction  drag  would  be  a  function  that 
predicts  the  percentage  of  attached  flow  and  laminar  flow.  Iterative  design  changes  could 
be  made  that  maximize  laminar  flow  and  minimize  separated  flow. 
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in.  CMARC  PANEL  CODE  THEORY 

A.  POTENTIAL  FLOW  PANEL  CODE  THEORY  (CMARC/PMARC) 

Potential  flow  theory  involves  the  superposition  of  sources  and  doublets  to 
generate  the  desired  flow  field  around  a  3D  body.  It  assumes  inviscid,  irrotational  and 
incompressible  flow.  As  such,  valid  solutions  are  only  obtained  at  low  Mach  numbers  and 
for  flow  fields  without  large  areas  of  separation. 

The  basic  concept  of  panel  methods,  as  outlined  by  Bertin  and  Smith  [Ref  3], 
requires  the  modeling  of  the  desired  3D  configuration  with  a  large  number  of  quadrilateral 
and  triangular  panels  representing  the  surface  of  the  aircraft.  A  series  of  sources,  doublets 
and  vortices  is  then  distributed  on  each  panel.  Superposition  allows  the  simultaneous 
computation  of  the  singularity  strengths  required  to  satisfy  flow  tangency  on  the  surface. 
The  inviscid,  irrotational  and  incompressible  flow  field  represented  by  the  superposition  of 
sources  and  doublets  satisfies  the  Laplace  equation: 


3.1 


Using  Green’s  Theorem,  the  potential  at  any  point  P  in  the  flow  is  represented  by: 

Where  (o - O.)  represents  the  potential  firom  the  doublet  distribution  and 
n  ■  (V<I>  -  VO,  )  represents  the  potential  from  the  source  distributions. 

CMARC  is  a  low  order  panel  code  that  assumes  constant  source  and  doublet 
strength  distributions  across  each  panel.  Figure  3.1  shows  a  panel  layout  for  a  generic  3D 
wing  fuselage  configuration.  It  is  important  to  note  that  for  a  3D  solution,  there  is  an 
equivalence  to  surface  doublet  and  surface  vortex  distributions.  CMARC  implements 
source  and  doublet  distributions. 
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Figure  3.1  Typical  Wing-Body  Panel  Code  Conflguration,  from  Ref.  (4]. 

As  mentioned  previously,  the  general  boundary  condition  imposed  is  tangential 
flow  at  the  surface.  CMARC,  as  outlined  in  Ref  [2],  allows  the  modification  of  this 
boundary  condition  on  individual  panels  or  groups  of  panels.  A  normal  surface  velocity 
distribution  may  be  specified  to  simulate  flow  into  or  out  of  ducts. 

In  order  to  produce  lifl;,  a  potential  flow  panel  code  requires  a  method  to 
implement  the  Kutta  condition.  As  noted  in  Anderson  [Ref  5],  the  Kutta  condition  at  the 
trailing  edge  implies  that  the  circulation,  F,  around  an  airfoil  is  such  that  the  flow  exits  the 
trailing  edge  smoothly.  In  addition,  the  velocities  leaving  the  top  and  bottom  surfaces  are 
finite  and  equal  in  magnitude  and  direction. 

Panel  codes  impose  the  Kutta  condition  by  the  shedding  of  wake  panels  along  the 
trailing  edges  or  separation  lines.  Wake  panels  are  similar  to  a  surface  panel  with  only  a 
doublet  distribution.  The  doublet  strength  of  the  attached  wake  panel  equals  the 
difference  in  doublet  strengths  of  the  two  adjacent  surface  panels. 


The  CMARC  core  panel  code  processing  engine  is  functionally  equivalent  to  the 
PMARC  panel  code  module.  The  implemented  equations  are  well  documented  by  Ashby 
et  al.  [Ref.  4].  The  PMARC  documentation  includes  a  wing-body  combination,  shown  in 
Figure  3.1,  evaluated  by  PMARC  with  good  correlation  to  experimental  data.  The  results 
are  shown  in  Figures  3.2  and  3.3.  In  addition,  Lambert  [Ref  6]  compared  PMARC  panel 
code  results  to  several  theoretical  and  experimental  test  cases  with  good  correlation  at  low 
angle-of-attack.  Sensitivity  to  wake  placement  is  highlighted  by  his  studies. 

Wake  positioning  can  have  a  large  influence  on  potential  flow  solutions.  A  wake  is 
obviously  attached  to  the  trailing  edge  of  wings  and  control  surfaces  with  sharp,  thin 
trailing  edges  to  produce  the  Kutta  condition.  However,  wake  positioning  on  streamlined 
fuselages,  missile  airframes  and  nacelles  is  more  of  an  art  than  science.  Recently,  Tuncer 
and  Platzer  [Ref  7]  investigated  generalized  wake  placement  techniques  for  cylindrical 
bodies  of  revolution  with  good  correlation  to  experimental  data  at  up  to  20  degrees  angle- 
of-attack.  The  techniques  are  used  in  this  study  with  success  for  the  verification  of 
CMARC  calculations  for  flow  over  an  inclined  6;  1  prolate  spheroid. 


Figure  3.2  Comparison  of  Experimental  Data  and  PMARC  Results  for  Two 
Spanwise  Stations  of  the  Wing/Body  (a  =  4°),  from  Ref.  [4J. 


Percent  fuselage  length 


o  exp.  (ref.  13) 
upper  surface 

•  exp.  (ref.  13) 
lower  surface 

PMj^C,  upper  surface 

•••  PMARC,  lower  surface 


Figure  3.3  Comparison  of  Experimental  Data  and  PMARC  Along  the  Fuselage 
Centerline  of  the  Wing/Body  Configuration  (a  =  4°),  from  Ref.  [4J. 


B.  CMARC  BOUNDARY  LAYER  ANALYSIS  THEORY 

CMARC  and  PMARC  use  the  same  two-dimensional  integral  method  to  calculate 
boundaiy  layer  characteristics  along  a  surface  streamline.  A  transition  model 
automatically  switches  from  laminar  to  turbulent  calculations.  The  developers  of  the 
PMARC  code  chose  a  2D  integral  routine  over  a  3D  finite  difference  grid  method  due  to 
speed  and  robustness  of  the  calculations  [Ref  4].  Building  a  finite  difference  grid  is  a 
difficult  and  time  consuming  process  requiring  the  user  to  develop  grids  over  complex  3D 
surfaces.  In  addition,  boundary  layer  calculation  times  can  easily  exceed  that  required  for 
the  basic  potential  flow  solution.  Reference  [8]  gives  a  good  outline  of  three-dimensional 
finite  difference  methods. 
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The  CMARC  and  PMARC  User’s  Guides  [Refs.  2  and  4]  contain  detailed 
discussions  on  the  development  of  the  CMARC/PMARC  boundary  layer  code  starting 
from  the  two-dimensional  momentum  equation: 


i2  +  H)* 


e  dU 

U  dr] 


3.3 


The  momentum  integral  equation  is  numerically  integrated  along  a  surface  streamline.  The 
derivation  leading  to  Equation  3.3  is  developed  in  Appendix  A  for  completeness. 

The  laminar  region  of  the  boundary  layer  is  modeled  by  numerically  integrating  the 
following  exact  differential  equation.  The  equation  is  solved  iteratively  through  numerical 
integration  along  a  streamline  starting  at  a  stagnation  point  [Ref  4]; 


+  2222g{K,fx)p{rjf  dr]  +  d{Q) 


^^(o)y 


3.4 


Where;  U  -  velocity  at  outer  edge  of  boundary  layer 

0  -  momentum  thickness 

V  dr] 

q  -  generalized  coordinate  along  a  streamline 

The  value  g(K,u)  is  based  on  exact  solutions  for  a  number  of  pressure 
distributions.  Initial  work  was  conducted  by  Thwaites  with  improvements  by  Curie 
[Ref  9]; 

j(^,//)  =  F.W-.uG.W-0.45  +  6X  3.5 


CMARC  uses  an  empirical  transition  model  based  on  the  average  pressure 
gradient,  K ,  for  predicting  laminar  to  turbulent  transition.  The  following  relations  are 
used  to  calculate  the  transition  point  [Ref  4]: 
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3.6 


\Kdn 


n-rii, 


Where  is  the  streamline  coordinate  at  instability.  And,  K  is  the  local  pressure  gradient 
at  boundary  layer  instability  [Ref.  4]: 


K  =  -0.4709  +  0.1 1066  *  ln(Reg)  +  0.0058591  *  In^(Re^) 
K  =  0.69412  -  0.23992  *  ln(Re^)  +  0.0205  *  In'(Re^) 


(0<Re5  <650) 
(650<Re^<10000) 


3.7 


The  local  Reynolds  number  at  transition  is  correlated  to  K  with  the  following  expressions 
[Ref  4]: 


K  =  -0.0925  +  0.00007  *  Re,  (o  <  Re,  <  750) 

K  =  -0.12571  +  0.000114286*Re,  (750  <  Re,  <  1  lOO)  3.8 

=  1.5938 1- 0.45543  *ln(Re,)  + 0.032534  *ln^(Re,)  (l  100  <  Re,  <  3000) 


At  transition,  the  initial  turbulent  shape  factor,  H,  is  given  by  the  following  empirical 
formula  that  is  a  fit  to  data  developed  by  Coles  [Ref  9]; 


H  = 


1.4754 
log, 0  (Re,) 


+  0.9698 


3.9 


Provisions  are  made  to  check  for  turbulent  reattachment  if  laminar  separation  is 
encountered.  At  laminar  separation,  a  point  calculation  is  made  to  determine  if  the 
boundary  layer  will  reattach.  If  reattachment  is  predicted,  the  boundary  layer  code 
immediately  switches  to  turbulent  calculations.  No  attempt  is  made  to  model  the  laminar 
separation  bubble  or  provide  a  transition  length.  After  laminar  separation  is  predicted,  the 
following  empirical  relations  are  used  to  determine  if  reattachment  occurs  [Ref  4]: 
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K  =  0.0227  -  0.007575  *  Re^  -  0.000001 1 57  *  Re/ 


(Re,  >125) 
(Re,  <  125) 


3.10 


K  =  -0.09 


The  boundary  layer  code  in  CMARC  uses  a  point  transition  model.  No  attempt  is 
made  to  model  a  more  representative  transition  length.  Turbulent  calculations  begin  at 
transition  using  the  Nash-Hicks  model  [Ref.  4].  Calculations  continue  along  the 
streamline  until  turbulent  separation  is  predicted  or  the  end  of  the  streamline  is  reached. 
No  boundary  layer  data  is  available  after  separation. 

The  authors  of  PMARC  caveat  that  their  boundary  layer  calculations  are  quite 
accurate  for  predominately  2D  flow  but  break  down  in  regions  of  large  cross  flow  near 
separation.  This  premise  will  be  first  tested  by  comparing  predominately  2D  flow  over  the 
inboard  region  of  a  high  aspect  ratio  wing  to  the  finite  difference  calculations  performed 
by  the  Naval  Postgraduate  School  Unsteady  Potential  Flow  Code  (UPOT).  Then  a 
comparison  is  made  to  experimental  data  for  flow  over  an  inclined  prolate  spheroid.  The 
6;  1  prolate  spheroid  is  chosen  because  of  the  availability  of  extensive  experimental  data. 
In  addition,  three-dimensional  flow  around  the  prolate  spheroid  is  similar  to  flow  around  a 
streamlined  slender  fuselage. 
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IV.  CMARC  VERIFICATION 


A.  VERIFICATION  OF  CMARC  AGAINST  PMARC 

The  first  step  in  CMARC  verification  is  comparison  with  NASA’s  PMARC  panel 
code.  CMARC  is  PMARC- 12  rewritten  in  C  from  FORTRAN.  Additionally,  CMARC  is 
compiled  for  hosting  on  an  IBM  compatible  PC.  Other  than  some  added  command  line 
functionality  and  significant  memory  management  improvements,  the  CMARC  basic  panel 
code  and  boundary  layer  routines  are  equivalent  to  PMARC  and  should  produce  the  same 
results.  However,  due  to  the  recent  fielding  of  CMARC,  the  author  felt  it  prudent  to  spot 
check  the  solutions  to  verify  equivalency. 

CMARC  and  PMARC  were  both  fed  an  identical  input  file  for  a  straight  NACA 
2415  wing  with  a  6.4  aspect  ratio  at  5  degrees  angle  of  attack.  The  input  file  is  listed  in 
full  in  Appendix  C.  Figures  4.1  and  4.2  show  CMARC  and  PMARC  pressure  coefficients 
cross-plotted  for  chordwise  and  spanwise  wing  stations  respectively.  The  results  overlay 
as  an  identical  match.  Figure  4.3  and  4.4  display  the  boundary  layer  calculations  for  skin 
fiiction  coefficient  and  displacement  thickness.  Again  the  results  overlay.  Integrated 
forces  and  moment  listings  were  also  identical.  From  this,  it  is  inferred  that  CMARC  and 
PMARC  produce  equivalent  results. 

Although  both  programs  produce  equivalent  results,  it  is  worthy  to  note  that  there 
are  occasionally  small,  insignificant  differences  in  floating  point  calculations  and  rounding. 
Some  results  differ  by  a  digit  in  the  sixth  decimal  place.  In  addition,  with  identical  input 
files,  there  can  be  a  difference  in  convergence  likelihood.  Occasionally,  PMARC  failed  to 
converge  when  CMARC  did.  Again,  floating  point  differences  are  the  most  likely  source 
of  the  disparity.  Regardless,  difference  in  the  rates  of  convergence  were  slight  and 
relatively  transparent  to  the  user.  However,  in  all  cases  CMARC  was  better  behaved  with 
a  higher  likelihood  of  convergence.  It  is  concluded  that  CMARC  and  PMARC  results  are 
interchangeable. 
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Figure  4.1  Comparison  of  CMARC  and  PMARC  Pressure  Coefficients  For 
Chordwise  Wing  Station. 
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Figure  4.2  Comparison  of  CMARC  and  PMARC  Pressure  Coefficients  For 
Spanwise  Wing  Station. 


Normalized  Chord  Position  -  x/c 

Figure  4.3  Comparison  of  CMARC  and  PMARC  Skin  Friction  Coefficient  for  an 
Upper  Wing  Surface  Streamline. 


Normalized  Chord  Position  -  x/c 

Figure  4.4  Comparison  of  CMARC  and  PMARC  Boundary  Layer  Displacement 
Thickness  for  an  Upper  Wing  Surface  Streamline. 
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B.  COMPARISON  OF  CMARC  AND  PMARC  PROCESSING  TIMES 


One  of  the  primary  metrics  in  determining  suitability  of  a  panel  code  hosted  on  an 
inexpensive  PC  is  processing  time.  CMARC’s  processing  speed  should  be  within  an  order 
of  magnitude  of  PMARC  hosted  on  the  NPS  Aeronautics  Department  Silicon  Graphics 
(SGI)  workstations  to  be  of  much  utility.  To  this  end,  processing  times  were  compared 
for  identical  input  models  ranging  from  200  to  1600  panels.  Comparisons  were  performed 
between  a  150  MHz/48  MB  Pentium  PC  and  two  configurations  of  networked  SGI 
Indigo^  workstations.  One  workstation  was  the  150  MHz/64  MB  Indigo^  (Viper)  running 
the  IRIX  5.3  operating  system  and  the  other  a  250  MHz/128  MB  SGI  Indigo^  (Aurora) 
workstation  running  IRIX  6.2.  With  both  workstations,  file  input/output  is  addressed 
through  the  network  to  the  server. 

Three  versions  of  PMARC  were  tested.  The  first  version,  “pmarcl2”  located  in 
the  local/usr^in,  was  compiled  in  FORTRAN  to  run  on  the  older  IRIX  5.3  operating 
system.  The  second  version,  “pmarc-inram,”  was  compiled  with  D)maniic  Linked 
Libraries  (DLLs)  to  run  on  the  new  IRIX  6.2  operating  system.  These  first  two  PMARC 
executable  codes  were  compiled  with  the  “in  RAM’  flag  selected  for  matrix  storage.  This 
considerably  reduced  hard  disk  accessing  time.  PMARC  either  keeps  all  the  matrices  in 
RAM  or  all  on  the  file  server  depending  on  whether  the  RAM  flag  is  when  compiled. 
Flexible  memory  management  allows  CMARC  to  fill  available  RAM  and  then 
automatically  spill  over  to  the  hard  drive.  This  reduces  user  memory  management 
requirements. 

A  third  version,  “pmarc  dll,”  was  compiled  with  DLLs,  but  the  matrix  storage  flag 
was  inadvertently  set  to  hard  drive  instead  of  RAM.  Processing  times  were  considerably 
longer  with  this  option  selected.  It  is  not  recommended  unless  the  computer  is  RAM 
limited. 

The  input  model  used  was  a  NACA  2415  finite  wing  with  four  time  steps.  The 
panel  density  was  varied  to  obtain  the  desired  panel  count.  Appendix  D  contains  a 
representative  input  file. 

The  processing  benchmarks  showed  that  CMARC,  hosted  on  a  PC,  is  significantly 
faster  than  all  three  versions  of  PMARC  hosted  on  the  networked  SGI  workstations. 
Table  4.1  summarizes  results  for  identical  models  ranging  from  200  -  1600  panels. 

Figure  4.5  is  a  plot  of  processing  time  vs.  panel  count  for  models  ranging  from 
200  to  1600  panels.  For  small  sized  models,  all  configurations  are  relatively  close  to  the 
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same  speed.  As  the  model  size  increases,  processing  time  increases  roughly  as  the  square 
of  model  size.  However,  as  model  size  increases,  all  versions  of  PMARC  on  the 
networked  SGI  workstations  become  significantly  slower  than  CMARC  on  the  PC.  This  is 
most  likely  due  to  the  slower  file  read/write  access  times  to  the  file  server.  The  version  of 
PMARC  with  the  matrix  storage  flag  set  to  hard  drive  required  considerably  more 
processing  time  than  the  two  versions  with  RAM  selected. 


Platform 

Pentium  PC 

SGI  Indigo^ 

SGI  Indigo^ 

CPU  /  RAM 

150  MHz/ 48  MB 

150  MHz /64  MB 

250  MHz  / 128  MB 

Program 

CMARC 

PMARC 

PMARC-DLL 

PMARC 

PMARC-DLL 

Panel  Count 

minisec 

minrsec 

minisec 

min:sec 

minisec 

200 

0:11 

0:12 

0:15 

N/A 

0:13 

400 

0:27 

0:43 

0:53 

N/A 

0:29 

800 

1:29 

2:46 

3:10 

N/A 

1:51 

1600 

5:54 

9:31 

17:30 

N/A 

10:04 

Table  4.1  CMARC  and  PMARC  Processing  Times  for  Models  Ranging  from  200 
to  1600  Panels. 

It  is  important  to  note  that  the  models  compared  in  this  study  only  differed  in  panel 
count.  Panel  count  is  not  the  only  factor  in  determining  processing  time.  The  number  of 
time  steps  selected,  solution  resolution,  convergence  rate  and  boundary  layer  calculations 
will  all  impact  processing  speed.  As  a  result,  the  times  presented  should  only  be  viewed  as 
representative  of  the  relative  impact  of  panel  density  and  not  as  the  time  required  to 
process  any  other  model  geometry. 

In  conclusion,  CMARC  hosted  on  a  dedicated  150  MHz  Pentium  PC  is 
significantly  faster  than  PMARC  hosted  on  a  similar  or  faster  networked  SGI  workstation. 
In  some  cases,  over  twice  as  fast.  Clearly,  executing  the  CMARC  panel  code  on  the  PC  is 
a  suitable  alternative  to  running  PMARC  on  the  SGI  workstations.  Low  cost  200-300 
MHz  Pentium  II  PCs  are  now  available  which  will  allow  further  reductions  in  CMARC 
processing  times. 
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Figure  4.5  Comparison  of  CMARC  and  PMARC  Processing  Times  for  Similar 
Finite  Wing  Models  Ranging  from  200  to  1600  Panels. 
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C.  COMPARISON  OF  CMARC  TO  THE  UPOT  BOUNDARY  LAYER  CODE 

As  a  first  step  in  investigating  CMARC  boundary  layer  calculations  and  utility, 
CMARC  results  are  compared  to  2D  calculations  from  the  NPS  Unsteady  Potential  Flow 
Code  (UPOT).  Although  the  potential  flow  solution  used  by  CMARC  for  the  boundary 
layer  calculations  is  strictly  a  3D  solution,  2D  flow  can  be  approximated  with  the  proper 
choice  of  geometry.  In  this  case,  flow  over  the  inboard  portion  of  a  high  aspect  ratio  (AR) 
wing  is  selected.  A  straight  NACA  2415  wing  with  AR=20  is  chosen  for  the  comparison. 
The  NACA  2415  is  the  same  section  used  in  the  FROG  UAV.  Boundary  layer  transition 
and  separation  points  are  compared  at  angles-of-attack  ranging  fi'om  0  to  20  degrees.  In 
addition,  boundary  layer  solution  sensitivity  is  investigated  over  four  Reynolds  numbers 
ranging  from  5.0x10^  to  6.0x10^. 

1.  UPOT  Boundary  Layer  Calculations 

The  NPS  UPOT  panel  code  was  developed  as  a  tool  to  assist  in  unsteady  flow 
visualization  over  two-dimensional  airfoils.  It  features  an  excellent  interactive  graphical 
user  interface  and  rapid  modeling  capabilities  [Ref  10].  Unlike  the  integral  momentum 
equations  used  by  CMARC  and  PMARC,  UPOT  implements  the  Cebeci-Keller  finite 
difference  boundary  layer  code.  The  algorithm  is  documented  by  Nowak  [Ref  11].  The 
UPOT  code  has  been  compared  to  experimental  data  for  a  range  of  airfoils  with  favorable 
results.  As  such,  it  is  considered  to  be  acceptable  to  benchmark  CMARC  results. 

2.  High  AR  Wing  Model 

A  high  aspect  ratio  NACA  2415  wing  is  modeled  to  evaluate  the  boundary  layer 
over  the  inboard  section  to  approximate  2D  flow.  CMARC’s  built-in  modeling  capability 
was  used  to  generate  a  finite  wing  with  dimensions  of  20  ft  wingspan  (b)  and  unit  chord 
(c)  yielding  an  aspect  ratio  of  20.  Fifty  chordwise  panels  are  distributed  over  the  top  and 
bottom  surface  in  a  full  cosine  distribution  and  10  panel  sections  in  a  spanwise  direction 
with  half  cosine  distribution.  There  are  600  panels  total,  including  the  enclosed  wing  tip, 
over  the  semi-span.  Figure  4.6  displays  a  semi-planform  view  of  this  configuration. 
Streamlines  are  placed  on  the  upper  and  lower  surfaces  of  the  inboard  root  panels.  The 
root  area  is  chosen  as  the  area  where  the  flow  is  nearly  two-dimensional  flow. 
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Upper  Surface 
Streamline 


Figure  4.6  Semi-Span  of  Finite  Wing  for  the  Approximation  of  Two-Dimensional 
Flow  Near  the  Root  (AR=20).  50  Chordwise  x  10  Spanwise  Panels. 

3.  .  Boundary  Layer  Results  and  Analysis  (CMARC  vs.  UPOT) 

CMARC  and  UPOT  boundary  layer  calculations  are  compared  for  the  FROG 
UAV  NACA  2415  airfoil.  Two  angles-of-attack  were  chosen  for  comparison.  The  first, 
-2°  or  zero  lift,  is  used  to  compare  transition  models.  The  second  angle-of-attack,  10°  is 
selected  for  comparison  to  the  10°  incidence  of  the  inclined  spheroid  discussed  in  a  later 
section.  A  comparison  for  Reynolds  numbers  ranging  fi-om  0.5x10^  to  6.0x10®  is  also 
performed  at  10°  to  investigate  boundary  layer  calculation  sensitivity  to  Reynolds  number. 
In  addition,  boundary  layer  transition  and  separation  locations  are  compared  at  angles-of- 
attack  ranging  from  0°  to  20°  at  Re=l  .0x10®. 

a.  Boundary  Layer  Transition 

The  shortcomings  of  the  point  boundary  layer  transition  model  coded  in 
CMARC  is  evident  when  compared  to  the  more  sophisticated  transition  length  model 
implemented  in  UPOT.  UPOT  uses  the  Michel  transition  onset  and  the  Chen-Thyson 
transition  length  models  [Ref.  11].  Figures  4.7  through  4.14  display  skin  friction 
coefficient  as  a  function  of  chordwise  location  (x/c)  for  the  upper  and  lower  surfaces  of  a 
NACA  2415  airfoil.  Results  for  four  Reynolds  numbers  ranging  from  0.5x10®  to  6.0x10® 
are  plotted  at  zero  lift  (-2°)  and  10°  angle-of-attack.  Boundary  layer  transition  will  be 
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discussed  first,  followed  by  boundary  layer  separation.  Finally,  differences  in  modeling  at 
the  stagnation  point  will  be  discussed. 

In  almost  all  cases,  CMARC  predicts  an  early  transition.  The  transition 
from  laminar  to  turbulent  boundary  layer  occurs  in  CMARC  as  a  sudden  jump  or  point 
transition.  The  UPOT  transition  length  model  provides  for  a  more  realistic  representation 
of  the  boundary  layer  physics.  Combined,  early  and  point  transition  result  in  higher  total 
skin  fiiction  drag  predicted  by  CMARC.  The  difference  is  most  pronoimced  at  the  lower 
Reynolds  numbers  associated  with  the  FROG  UAV.  At  Re=0.5xl0^  and  zero  lift, 
CMARC  overpredicts  skin  fiiction  drag  by  approximately  40%  on  the  upper  surface  and 
20%  on  the  lower  surface.  Although  skin  friction  drag  may  be  a  relatively  small  portion  of 
the  total  drag,  airframe  manufacturers  go  to  great  lengths  to  refine  models  to  accurately 
predict  it.  A  few  percentage  points  of  error  can  cause  the  aircraft  to  meet  or  miss 
performance  goals. 

Despite  the  differences  in  transition  modeling,  CMARC  accurately  predicts 
the  skin  fiiction  coefficient  with  respect  to  UPOT.  When  comparing  laminar  to  laminar 
and  turbulent  to  turbulent  regions  in  Figures  4.7  through  4.10  (zero  lift  plots),  the  skin 
fiiction  coefficients  are  a  close  match.  This  indicates  that  an  adjustment  in  the  CMARC 
model  delaying  transition  could  provide  more  accurate  results. 

As  another  comparison  of  boundary  layer  calculations,  displacement 
thickness  (6*)  is  displayed  in  Figures  4.15  through  4.22  as  a  function  of  chord  position 
(x/c)  for  zero  lift  (-2°)  and  10°  angle-of-attack.  In  general,  CMARC  and  UPOT  predict 
similar  trends  in  5*.  The  final  displacement  thickness  is  a  good  relative  indication  of  total 
skin  fiiction  drag.  CMARC  always  predicts  a  greater  6*  and  thus  more  drag.  This  is  in 
keeping  with  the  previous  observations  indicating  higher  integrated  skin  fiiction  forces. 

b.  Separation 

Boundary  layer  separation  is  indicated  in  Figures  4.7  through  4.14  by  a 
zero  or  negative  skin  fiiction  coefficient.  In  all  cases,  CMARC  slightly  overpredicts  the 
extent  of  attached  flow.  Again,  the  differences  are  most  significant  at  the  lower  Reynolds 
numbers. 

Figures  4.23  and  4.24  display  transition  and  separation  points  for  the 
NACA  2415  as  a  function  of  angles-of-attack  ranging  from  0°  to  20°.  The  data  is  for 
Re=  1.0x10^  which  is  close  to  the  FROG  UAV  high  speed  cruise  at  Re=929,000.  On  both 
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Figure  4.7  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  zero  lift  (a=-2°)  and  Re=0.5xl0^ 


Figure  4.8  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  zero  lift  (a=-2°)  and  Re=1.0xl0^ 
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Figure  4.9  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  zero  lift  (a— 2°)  and  Re=3.0xl0^. 


Chord  Position  -  x/c 

Figure  4.10  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  zero  lift  (a— 2°)  and  Re=6.0xl0®. 
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Figure  4.11  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  a=10°  and  Re=0.5xl0^ 


Chord  Position  -  x/c 

Figure  4.12  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  a=10°  and  Re=1.0xl0^ 
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Figure  4.13  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  a=10°  and  Re=3.0xl0®. 


Figure  4.14  Comparison  of  CMARC  and  UPOT  Skin  Friction  Coefficient  (Cf)  for 
NACA  2415  at  a=10'’  and  Re=6.0xl0^ 
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Figure  4.15  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  zero  lift  (a=-2°)  and  Re=0.5xl0*. 


Chord  Position  -  x/c 

Figure  4.16  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  zero  lift  (a— 2°)  and  Re=1.0xl0^ 
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Figure  4.17  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (8*)  for  NACA  2415  at  zero  lift  (a=-2°)  and  Re=3.0xl0^. 


Chord  Position  -  x/c 

Figure  4.18  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  zero  lift  (a=-2°)  and  Re=6.0xl0^ 
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Figure  4*19  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  a=10°  and  Re=0*5xl0^ 


Figure  4.20  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  a=10°  and  Re=1.0xl0^ 
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Figure  4.21  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  a=10°  and  Re=3.0xl0®. 
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Figure  4.22  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Displacement 
Thickness  (5*)  for  NACA  2415  at  a=10°  and  Re=6.0xl0®. 
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i.23  Comparison  of  CMARC  and  UPOT  Boundaiy  Layer  Transition  and 
Separation  Points  for  the  Upper  Surface  of  a  NACA  2415  Airfoil  at 
Re=1.0xl0^  from  0°  to  20°  AO  A. 
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Figure  4.24  Comparison  of  CMARC  and  UPOT  Boundary  Layer  Transition  and 
Separation  Points  for  the  Lower  Surface  of  a  NACA  2415  Airfoil  at 
Re=1.0xl0®  from  0°  to  20°  AOA. 
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the  upper  and  lower  surface,  CMARC  clearly  provides  correct  trends  for  both  the 
transition  and  separation  points.  However,  as  seen  at  zero  lift  in  Figures  4.7  through  4.10, 
CMARC  always  predicts  an  early  transition  and  late  separation. 

Despite  the  inaccuracies  in  transition  and  separation  points,  CMARC 
boundary  layer  calculations  remain  useful.  A  low  order  panel  code  is  unlikely  going  to  be 
used  for  performance  calculations.  Instead,  it  is  more  useful  as  a  design  tool.  It  allows  for 
rapidly  visualizing  the  trend  in  transition  and  separation  points  with  minor  changes  in 
configuration. 

A  word  of  caution  is  advised  when  total  skin  friction  drag  is  integrated.  A 
design  change  could  be  implemented  that  reduces  overall  skin  fnction  drag  but  neglects 
large  increases  in  pressure  drag.  In  other  words,  one  could  reduce  skin  friction  drag,  but 
fail  to  realize  earlier  separation  is  taking  place.  The  net  result  is  a  small  reduction  in  skin 
fiiction  drag  that  is  more  than  offset  by  a  large  increase  in  separation  pressure  drag. 
Extending  the  extent  of  attached  flow  should  always  be  considered  preferable  to  reducing 
overall  integrated  skin  fnction  drag. 

a  Skin  Friction  Confident  near  the  Stagnation  Point 

Another  major  difference  between  the  integral  boundary  layer  code  in 
CMARC  and  the  finite  difference  code  in  UPOT  is  highlighted  at  the  stagnation  point.  In 
Figure  4. 1 1,  both  codes  locate  the  stagnation  point  on  the  lower  surface  at  x/c=0.025  for 
the  NACA  2415  airfoil  at  10  degrees  angle-of-attack.  However,  it  is  clear  that  the 
CMARC  skin  fnction  coefficient  starts  at  0.0002,  a  small  number  approaching  zero 
asymptotically,  while  the  UPOT  skin  fnction  shoots  out  of  the  top  of  the  chart  in  excess 
0.7,  a  relatively  large  number  approaching  +oo  asymptotically. 

From  boundary  layer  theory,  it  is  know  that  the  skin  fiiction  coefficient  is 
inversely  proportional  to  the  square  root  of  the  local  Reynolds  number  or: 


C. 


4.1 


At  the  stagnation  point,  Cf  approaches  +<x).  The  finite  difference  code  in 
UPOT  correctly  models  this  trend.  On  the  other  hand,  CMARC  implements  a  discrete 
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integration  of  the  following  exact  differential  laminar  skin  friction  calculation: 


Where;  U  -  velocity  at  outer  edge  of  boundary  layer 

0  -  momentum  thickness 
g  -  empirical  parameter 

V  dr} 

T|  -  generalized  coordinate  along  a  streamline 


At  T|  -0,  the  momentum  thickness  starts  at  zero  and  builds  rapidly  from  the 
stagnation  point.  Thus,  the  momentum  integral  equation  reduces  to; 


dr] 


at  the  stagnation  point. 


4.3 


The  integral  solution  for  Q  starts  at  zero  and  rises  quickly  until  the  integral  portion  of 
Equation  4.2  dominates. 

The  incorrect  modeling  of  Cf  near  the  stagnation  point  is  considered  minor 
due  to  its  local  nature  at  the  stagnation  point.  When  skin  friction  is  integrated  over  the 
entire  surface  the  differences  are  bound  to  be  relatively  small.  In  addition,  when  integrated 
into  a  force,  the  errors  in  Cf  at  the  stagnation  point  tend  to  cancel  out.  Close  to  the 
stagnation  point  on  either  side,  skin  friction  forces  are  opposite  in  direction. 

D.  COMPARISON  OF  CMARC  TO  INCLINED  PROLATE  SPHEROID 
EXPERIMENTAL  DATA 

In  the  previous  section,  model  geometry  was  selected  to  produce  predominantly 
two-dimensional  flow.  In  this  section,  CMARC  pressure  distributions  and  integral 
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boundary  layer  data  are  investigated  for  a  model  geometry  that  produces  largely  three- 
dimensional  flow.  For  comparison,  a  suitable  experimental  test  case  was  found  in 
AGARD  AR-303:  A  Selection  of  Experimental  Test  Cases  for  the  Validation  of  CFD 
Codes  [Ref  12],  Case  number  C-2,  entitled  “Three-Dimensional  Boundary  Layer  and 
Flow  Field  Data  of  an  Inclined  Prolate  Spheroid”  was  selected.  A  6;  1  prolate  spheroid 
approximates  a  typical  streamlined  fuselage.  The  data  set  was  ordered  from  AGARD 
through  the  NASA  Center  for  Aerospace  Information  (CASI). 

A  complete  data  set  for  all  test  cases  in  AGARD  AR-303  was  available  for  a 
nominal  charge  of  $59.00  through  NASA’s  CASI  publications  office.  Ordering 
information  inside  the  rear  cover  of  the  publication  proved  to  be  accurate  and  useful.  The 
data  arrived  on  nine  PC  formatted  high  density  disks.  After  copying  the  desired  data  sets 
to  the  hard  drive,  each  file  is  self  extracting  through  a  built-in  DOS  decompression 
program.  Detailed  instructions  are  printed  in  the  back  section  of  AR-303. 

1.  Inclined  6:1  Prolate  Spheroid  -  AGARD  AR-303  Case  C-2 

AGARD  AR-303  test  case  number  C-2  contmns  pressure  coefficient  and  skin 
friction  distributions  for  a  6: 1  prolate  spheroid  inclined  to  the  flow  field.  Table  4.2  lists 
the  test  conditions  for  which  data  are  available. 

Test  case  I  was  chosen  for  comparison  to  CMARC  output.  At  10°  angle-of- 
attack,  some  separated  flow  was  expected  which  would  provide  a  good  comparison  to 
CMARC  integral  boundary  layer  separation  points.  The  only  drawback  to  this  test  case  is 
the  forced  transition  at  X/2a  =  0.20.  Natural  transition  would  have  been  more  desirable 
for  comparison  to  the  CMARC  transition  model.  The  test  cases  at  30°  angle-of-attack  are 
deemed  to  have  too  much  separated  flow  to  provide  a  meaningful  comparison  to  a 
CMARC  potential  flow  solution. 

a.  Wind  Tunnel  Experimental  Set-up 

Figure  4.25  contains  a  diagram  of  the  experimental  set-up  for  the  6:1 
prolate  spheroid  performed  by  Kreplin  in  the  DLR  Gottingen  three  meter  Low  Speed 
Wind  Tunnel  (NWG).  Of  note,  the  wind  tunnel  test  section  is  of  the  Gbttingen  type  with 
closed  return  and  open  test  section.  No  corrections  are  applied  to  the  data. 
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PARAMETER 

CASE  I 

CASE  II 

CASE  III 

Mach  Number 

0.16 

0.13 

0.23 

Reynolds  Number 

7.7  X  10"* 

6.5  X  10® 

43.0  X  10® 

Incidence 

10.0° 

29.7° 

30.0° 

Transition 

tripped  at  X/2a  =0.20 

free 

free 

Table  4.2  AGARD  AR-303  Test  Conditions,  from  Ref.  [12]. 


Figure  4.26  shows  the  configuration  for  the  6:1  prolate  spheroid  wind 
tunnel  model.  The  2.4  meter  long  model  contains  42  pressure  taps  located  along  an  axial 
meridian.  The  model  can  be  rotated  axially  in  50  steps  through  just  over  180  degrees 
providing  in  excess  of 2000  pressure  readings  over  half  the  surface.  With  yaw  angle  set  to 
zero,  symmetry  is  assumed  for  the  other  half.  In  addition  to  pressure  ports,  surface  hot 
film  sensors  are  located  at  12  axial  positions  for  the  measurement  of  wall  shear  stress. 
Wall  shear  stress  is  normalized  by  dynamic  pressure  to  provide  skin  friction  coefficient 
(Cf).  Once  again,  the  measurements  are  provided  for  approximately  50  rotation  angles, 
providing  coverage  of  half  the  surface. 

In  addition  to  pressure  and  skin  friction  coefficients,  boundary  layer 
velocity  profiles  and  flow  field  mean  velocity  vectors  are  available  at  several  axial 
locations.  Although  not  used  in  this  investigation,  the  data  would  prove  useful  for  more 
detailed  studies. 

Unfortunately,  the  wind  tunnel  set-up  was  not  instrumented  for  loads.  As 
will  be  discussed  in  the  next  section,  the  number  of  pressure  and  skin  friction 
measurements  was  deemed  to  be  sufficient  to  allow  the  integration  of  local  forces  to 
provide  a  reasonable  calculation  of  lift,  drag  and  pitching  moments. 

b.  Experimental  Data 

Two  data  files  from  test  case  C-2  at  a=10.0°  are  used  for  comparison  to 
CMARC  data.  The  first  file,  “cpl0nwg.dat,”  contains  pressure  coefficient  listed  as  a 
function  of  axial  location  (X/2a)  and  circumferential  angle  (<()).  For  each  circumferential 
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Figure  4.25  Inclined  6:1  Prolate  Spheroid  Model  in  the  DLR  Gottingen  Three 
Meter  Low  Speed  Wind  Tunnel  (NWG) ,  from  Ref.  [12]. 


Figure  4.26  Prolate  Spheroid  Wind  Tunnel  Model  Configuration,  from  Ref.  [12]. 
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angle  all  the  successive  axial  location  pressure  coefficients  were  listed.  It  is  more  common 
to  plot  pressure  distribution  as  a  function  of  circumferential  angle  at  a  given  axial  station. 
The  data  file  was  rearranged  using  the  MATLAB  M-file  listed  in  Appendix  E. 

The  second  data  file,  “cflOnwg.dat,”  contains  skin  fiiction  listed  as  a 
function  of  axial  location  (X/2a)  and  circumferential  angle  (<j)).  The  data,  listed  in  two 
columns,  were  reordered  to  one  column  for  ease  of  plotting. 

c.  Integration  of  Local  Forces  to  Provide  Lift  Drag  and  Pitching 

Moment 


The  experimental  set-up  did  not  include  balance  measurement  of  forces. 
However,  it  was  deemed  that  the  2000+  pressure  and  500+  skin  fiiction  measurements 
would  be  sufficient  to  allow  the  integration  of  measurements  over  the  surface  of  the 
prolate  spheroid  for  an  approximation  of  total  force  and  moment  coefficients.  The 
following  equations  were  developed  to  provide  integrated  pressure  and  fiiction  force 
coefficients.  Symmetry  is  assumed.  Appendix  B  outlines  the  development  of  these 
relations.  Appendix  C  lists  the  MATLAB  program  which  implements  the  discrete 
integration. 

The  pressure  force  coefficients  normalized  by  S-jdf  and  c  =  2b  are 
yielded  by  discretely  integrating  the  following  equations  in  a  cylindrical  coordinate  system: 
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{qSp^) ■  Axz  / 2a 

)=1  y=i 

m  n 

Mp  =  2^ [(^coCp/if) \xj2a-k-zJ2a- F)]a^^ Ax,.  / 2a 

i=l  M  ■' 


4.3 

4.4 

4.5 


Where  the  surface  unit  normal  is  given  by: 

m  ^  sin(^)  _ 


Unit  Normal:  n 


+1  + 1  +1 


_co^ 

J  ~  I — I - 


4.6 
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The  skin  fnction  coefficients  normalized  by  S  =  and  c  =  2h  are 
yielded  by  discretely  integrating  the  following  equations  in  a  cylindrical  coordinate  system; 


r 

c. 


,  ^SF  =  2^  S  /  la 


;=1  y=l 

m  n 


^SF  =  /  A^^  AX;  /  2a 


4.7 

4.8 


=  % ,  Msf=2±±  \[qSfry)  \-xJ2a-k+zJ  2a.  / )  /  2a  4.9 

1=1  y=i  *■ 


1=1  y=i 

m  /I 


Where  the  unit  surface  velocity  vector  is  given  by: 


cos(r)  _ 
v=  ,  + 

^|m^  + 1 


wsmi 


(<zi)cos(r) 


+  1 


+  cosi 


Wsin(r) 


J+ 


/«cos(^)cos(x)  .  /  \ 

+  sin(^jsin(;^j 


+1 
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The  surface  and  local  slope  of  a  prolate  spheroid  comes  from  the  following  relations: 


Prolate  Spheroid;  ^  ^  ^ 


dr 

slope  /w  =  —  = 
dx 


bx 


a\l- 


a‘ 


4.11 


Note:  The  forces  are  summed  over  half  the  spheroid,  ^  =  0  — >  1 80° ,  and  doubled. 

The  y-direction  forces  and  the  roll  and  yaw  moments  are  neglected  zero  due  to  symmetry. 


2.  CMARC  Model  of  6:1  Prolate  Spheroid 

A  40x20  panel  model  of  the  2.4  meter  6: 1  prolate  spheroid  wind  tunnel  model  was 
created  with  LOFTSMAN.  The  right  half  surface  was  modeled  with  symmetry  around  the 
y=0  plane.  Appendix  F  contains  a  printout  of  the  LOFTSMAN  input  file  which  includes  a 
fore/aft  wake.  The  patch  was  created  with  40  axial  and  20  semi-circumferential  panels. 
Full  cosine  compression  was  used  to  bunch  panels  at  the  leading  and  trailing  edge. 
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After  creating  the  patch  in  LOFTSMAN,  it  was  decided  that  a  doubling  of 
circumferential  panel  count  would  increase  wake  placement  flexibility.  The  CMARC  input 
file  was  modified  to  create  40  circumferential  panels  by  setting  TNPC=40  in  the  break 
point  input  field  for  each  cross  section. 

Figure  4.27  is  a  POSTMARC  rendering  of  the  final  1600  (40x40)  panel 
configuration.  The  input  file  takes  advantage  of  the  plane  of  symmetry  capability  built  into 
CMARC.  It  calculates  just  half  a  model  symmetric  around  the  y=0  plane  of  symmetry 
provided  there  is  zero  side  slip. 


Figure  4.27  CMARC  40x40  Prolate  Spheroid  Model  Rendered  with  POSTMARC. 


3.  Data  Extraction 

Pressure  coefficient  data  are  extracted  using  the  “postprolate.exe”  FORTRAN  file 
listed  in  Appendix  G.  This  program  extracts  data  from  a  CMARC  or  PMARC  output  file 
(DATA6)  for  a  range  of  panel  numbers  and  places  them  in  a  separate  plot  input  file. 
CMARC  output  files  are  transferred  to  the  SGI  workstations  for  data  extraction  using  the 
Windows  3.1  FTP  program.  Results  are  then  plotted  against  experimental  data  with  any 
x-y  plotting  program  (xmgr). 
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4. 


Prolate  Spheroid  Pressure  Distribution 


CMARC  and  experimental  pressure  coefficients  are  compared  at  10  degrees  angle- 
of-attack  and  Re=7.7xl0^.  Results  are  displayed  as  a  function  of  axial  station,  x/2a,  and 
circumferential  angle,  (j),  in  Figure  4.28.  Circumferential  angle  is  measured  starting  from 
the  lower  centerline  of  the  model.  CMARC  generated  potential  flow  pressure  coefficients 
over  the  forward  60%  of  the  prolate  spheroid  closely  match  experimental  results.  Of  note, 
there  is  a  constant  bias  between  the  two  sets  of  data. 

The  divergence  between  CMARC  and  experimental  results  aft  of  x/2a=0.60 
indicates  flow  separation  over  the  top  portion  of  the  prolate  spheroid.  It  is  clear  that  a 
potential  flow  solution  without  wakes  does  a  poor  job  of  predicting  pressure  distribution 
over  regions  with  separated  flow. 

To  model  the  flow  separation,  wakes  were  added  to  the  CMARC  model.  Tuncer 
and  Platzer’s  research  [Ref  7]  indicates  that  proper  wake  placement  can  produce  a  close 
match  between  panel  code  and  experimental  results  for  slender  bodies  of  revolution  for 
angles-of-attack  up  to  20  degrees.  They  concluded  that  a  circumferential  wake  placement 
angle  of  144  degrees  on  an  ogive  cylinder  body  provides  the  closest  match  for  force  and 
moment  coefficients. 

A  series  of  wakes  were  placed  at  several  circumferential  angles  ranging  from  117 
to  162  degrees.  The  wakes  run  fore-aft  from  x/2a=0.50  to  a  wake  separation  ring  at 
x/2a=0.99.  Results  are  plotted  in  Figure  4.29.  A  wake  angle  of  1 17  degrees  produced  the 
closest  average  match  to  the  experimental  results.  Figure  4.30  shows  the  final  wake 
configuration. 
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Pressure  Coefficient  -  Cp  Pressure  Coefficient  -  C  Pressure  Coefficient  - 


Circumferential  Angle  -  ^  (deg)  Circumferential  Angle  -  ^  (deg) 


Circumferential  Angle  -  ()>  (deg)  Circumferential  Angle  -  <|)  (deg) 


Circumferential  Angle  -  <ji  (deg)  Circumferential  Angle  -  ((i  (deg) 


Figure  4.28  CMARC  Potential  Flow  (No  Wakes)  Pressure  Distribution 
Compared  to  Experimental  Data,  after  Ref.  [12]. 
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Pressure  Coefficient,  Cp  Pressure  Coefficient,  C 
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Figure  4.29  CMARC  Pressure  Distribution  with  Wake  Angles  Ranging  from  117 
to  162°  Compared  to  Experimental  Data,  after  Ref.  [12]. 


Figure  4.30  POSTMARC  Views  of  CMARC  Model  with  117°  Wake  Separation 
Line  Running  Aft  from  x/2a=0.5  to  x/2a=0.99. 
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Coefficients  for  normal  (Cn),  axial  (Ca),  lift  (Cl),  drag  (Cd),  and  pitching  moment 
(Cm)  are  compared  to  experimental  forces  in  Table  4.3  for  a  circumferential  wake  angle  of 
117  degrees.  CMARC  automatically  outputs  the  pressure  load  coefficients  in  both  wind 
and  body  axes.  Skin  fnction  forces  are  calculated  using  POSTMARC  and  will  be 
discussed  in  a  later  section.  The  experimental  results  are  from  integrated  pressure  forces 
using  the  method  outlined  in  Appendix  B.  The  coefficients  are  normalized  by  maximum 
diameter  and  cross  sectional  area.  A  wake  angle  of  1 17  degrees  produces  a  close  match 
to  experimental  results  for  Cn,  Cl  and  Cm-  As  expected,  the  axial  and  drag  coefficients  are 
off  considerably  from  experimental  data. 


Force  Origin 

Force/Moment 

Coefficient 

Experimental 

AGARD  303-Kreplin 

CMARC 

e,sii7* 

%  Difference 
(CMARC-exp)/exp 

Cn 

0.1924 

0.1816 

-5.6% 

Pressure 

Ca 

0.0026 

0.0411 

1480.8% 

Forces 

Cl 

0.1890 

0.1717 

-9.2% 

Cd 

0.0359 

0.0720 

100.6% 

Cm 

0.9009 

0.9003 

-0.1% 

Table  4.3  Comparison  of  Integrated  Experimental  Pressure  Forces  to  the  CMARC 
Model  with  117®  Wake  Placement  Angle,  after  Ref.  [12]. 

It  is  concluded  that  a  pure  potential  flow  solution  over  a  streamlined  body  at  10° 
angle-off-attack  will  fail  to  predict  substantial  regions  of  flow  separation.  However, 
pressure  distributions  over  bodies  with  substantial  flow  separation  can  be  approximated  by 
proper  wake  distribution.  As  outlined  by  Tuncer  and  Platzer  [Ref.  7],  a  wake  separation 
angle  of  144°  is  a  good  starting  point. 

5.  Boundary  Layer  Separation  Locations 

Next,  CMARC  boundary  layer  calculations  were  visualized  to  see  how  well 
CMARC  predicted  separation  points  for  the  inclined  prolate  spheroid.  As  reported  in  the 
section  on  the  NACA  2415  finite  wing,  predicted  boundary  layer  separation  points  from 
CMARC  matched  those  predicted  by  the  NPS  UPOT  code  fairly  well,  especially  at  higher 
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Reynolds  numbers.  In  this  case,  the  boundary  layer  points  are  compared  to  experimental 
data  at  the  same  10°  angle-of-attack  over  the  three-dimensional  prolate  spheroid. 

Sixty-six  streamlines  for  boundary  layer  calculations  were  placed  on  the  CMARC 
model  at  locations  corresponding  to  experimental  data  points.  Appendix  F  contains  the 
input  file.  CMARC  only  predicted  separation  over  the  very  aft  end  of  the  body.  A 
separation  point  is  best  visualized  with  POSTMARC  by  selecting  the  on-body  streamline 
boundary  layer  thickness  or  shape  factor  functions.  The  separation  point  is  indicated  at 
the  last  downstream  point  on  the  streamline.  It  is  important  to  note  that  if  one  visualizes 
streamline  pressure  coefficient,  velocity  or  Mach  number,  the  streamline  will  travel  all  the 
way  to  the  aft  stagnation  point.  In  other  words,  to  visualize  a  separation  point, 
phenomena  derived  from  the  boundary  layer  calculations  and  not  the  streamline 
calculations  must  be  selected  for  visualization. 

Figure  4.31  displays  the  streamline  separation  points  on  the  aft  end  of  the  prolate 
spheroid.  CMARC  boundary  layer  separation  points  are  compared  to  experimental  data 
as  a  function  of  axial  location  and  circumferential  angle  in  Figure  4.32.  It  is  to  be 
expected  that  the  2D  code  implemented  in  CMARC  fails  to  accurately  predict  separation 
regions  over  streamlined  bodies  of  revolution  with  large  cross  flow  velocities. 
Nevertheless,  these  results  help  to  quantify  the  differences. 


Figure  4.31  POSTMARC  Visualization  of  CMARC  Predicted  Separation  Points  on 


the  Aft  End  of  the  Prolate  Spheroid  Model  (No  Wakes). 
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Figure  4.32  Comparison  of  CMARC  Predicted  Separation  Line  to  Experimental 
Data,  after  Ref.  [12]. 

6.  Boundary  Layer  Skin  Friction  Coefllcient 

CMARC-computed  skin  friction  coefficients  were  compared  to  experimental  data 
obtained  from  hot  film  sensors  [Ref  12],  Sixty-six  streamlines  were  placed  through 
panels  on  the  CMARC  model  corresponding  to  skin  fiiction  data  pints.  Data  were 
extracted  manually  from  the  CMARC  output  filew.  Data  are  plotted  at  six  axial  locations 
as  a  function  of  circumferential  angle  in  Figure  4.33.  The  wind  tunnel  model  has  a 
transition  strip  located  at  x/2a=0.20.  All  CMARC  boundary  layer  calculations  are  based 
on  a  built-in  transition  model.  There  are  no  provisions  for  specifying  the  transition 
location  in  CMARC. 
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Skin  Friction  Coeff.  Skin  Friction  Coeff.  Skin  Friction  Coeff. 
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Figure  4.33  Comparison  of  CMARC  Predicted  Skin  Friction  Coefficient  to 
Experimental  Data,  after  Ref.  [12]. 
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For  the  two  axial  locations  in  front  of  the  transition  strip,  x/2a=0.05  and  0.14, 
CMARC  predicts  a  mix  of  laminar  and  turbulent  flow.  A  laminar  boundary  layer  is 
indicated  by  the  data  points  where  Cf  <  0.002.  Experimental  data  indicate  strictly  laminar 
flow  for  these  axial  locations.  CMARC  streamlines  passing  through  each  circumferential 
location  travel  a  unique  path  across  different  panel  geometry  from  the  forward  stagnation 
point  to  the  point  of  interest.  Being  an  integral  two-dimensional  boundary  layer  method, 
CMARC’s  empirical  transition  formula  predicts  separation  for  some  of  the  streamlines  and 
laminar  flow  for  the  others.  In  general,  CMARC  over-predicts  skin  friction  drag  in  this 
region  due  to  the  mixed  flow.  If  CMARC  correctly  predicted  all  laminar  flow,  the  results 
would  be  close  to  experimental  results. 

Aft  of  the  transition  strip  at  x/2a=0.20,  experimental  data  indicate  fully  turbulent 
flow  as  expected.  CMARC  predicts  turbulent  flow  for  all  but  the  lower  streamline  which 
has  a  low  adverse  pressure  gradient.  At  x/2a=0.3 1  and  0.48,  computed  skin  fiiction  is 
accurate  to  within  25%.  Aft  of  x/2a=0.48,  CMARC  results  are  less  meaningful  due  to  the 
large  region  of  separated  flow. 

7.  Integrated  Skin  Friction  Forces 

POSTMARC  version  1.17.3  contains  functionality  for  performing  integrated  skin 
friction  calculations.  When  a  CMARC  model  is  processed  with  the  “-p”  command  line 
switch,  a  file  with  a  “.pm”  extension  is  created  with  the  information  necessary  for 
POSTMARC  to  perform  boundary  layer  calculations.  POSTMARC  then  places 
streamlines  on  every  panel,  performs  boundary  layer  calculations  and  integrates  the  skin 
fiiction  loads.  Experimental  data  is  integrated  as  outlined  in  Appendix  B. 

Integrated  skin  fiiction  forces  for  the  prolate  spheroid  model  without  wakes  are 
compared  to  experimental  data  in  Table  4.4.  Normal,  axial,  drag  and  pitching  moment 
coefficients  were  all  within  40%  of  the  rough  estimate  provided  by  integrating  the 
experimental  data.  This  is  in  keeping  with  the  observations  from  Figure  4.33.  The  lift 
coefficient  produced  due  to  skin  fiiction  is  so  small  that  comparisons  between 
experimental  and  CMARC  data  are  meaningless. 
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Force  Origin 

Force/Wloment 

Coefficient 

Experimental 

AGARD  3034<replin 

CMARC 

e«=ii7° 

%  Difference 

(CMARC-«xp)/exp 

Cn 

0.0102 

0.0071 

-30.6% 

Skin  Friction 

Ca 

0.0610 

0.0376 

-38.4% 

Forces 

Cl 

-0.0006 

0.0004 

-166.7% 

Cd 

0.0618 

0.0376 

-39.2% 

Cm 

0.0022 

0.0019 

-12.4% 

Table  4.4  Comparison  of  Integrated  Experimental  Skin  Friction  Forces  to  the 
CMARC  Model  without  Wakes,  after  Ref.  [13]. 

8.  Total  Integrated  Forces 

As  a  final  comparison  of  CMARC  results  to  experimental  data,  the  summed 
pressure  and  skin  fiiction  force  coefficients  are  presented  in  Table  4.5.  A  simple  fore/aft 
wake  running  from  x/2a=0.5  to  a  partial  ring  wake  at  x/2a=0.99  provides  good  results  for 
all  but  the  axial  and  drag  coefficients.  It  is  concluded  that  CMARC,  with  proper  wake 
selection,  will  provide  meaningful  force  and  moment  coefficients  for  the  development  of 
stability  derivative  data.  Results  for  drag  coefficient  are  less  meaningful  and  should  be 
avoided  for  performance  calculations. 
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V.  AERODYNAMIC  MODEL  OF  THE  FROG  UAV 


A.  BACKGROUND 

The  Naval  Postgraduate  School  Aeronautics  Department  is  integrating  UAV 
hardware  and  software  to  demonstrate  autonomous  flight,  trajectory  tracking  and 
automatic  landing.  A  core  requirement  for  flight  control  law  development  is  a  valid 
aerodynamic  truth  model  for  the  UAV  airframe.  A  panel  code  model  of  the  FROG  UAV 
is  one  method  for  estimating  many  of  the  stability  derivatives  required  for  an  aerodynamic 
truth  model.  This  development  effort  concentrates  on  finding  the  Clo  and  Cma  longitudinal 
stability  derivatives  followed  by  the  Cyp,  Cip  and  Cnp  lateral-directional  stability 
derivatives.  A  future  study  -will  continue  the  development  for  rate  damping  and  control 
effectiveness  derivatives. 

Panel  code  modeling  utility  goes  beyond  the  development  of  aerodynamic 
coefficients.  Flight  control  systems  require  accurate  pitot-static  and  angle-of-attack 
sensor  inputs.  CMARC  accurately  solves  on-body  static  pressure  distributions  and  oflf- 
body  flow  velocities  over  the  predominately  attached  flow  fields  of  fuselage  fore  bodies. 
In  this  study,  correction  curves  are  generated  for  static-pressure  source  and  angle-of- 
attack  probe  position  errors. 


B.  FROG  UAV  DESCRIPTION 

The  FROG  UAV  is  a  small  single  engine  flight  test  vehicle  used  for  autonomous 
flight  research  by  the  Naval  Postgraduate  School  Aeronautics  Department.  The  aircraft 
was  originally  designated  the  FOG-R  by  the  U.  S.  Army.  It  was  designed  as  a  small 
lightweight,  battlefield  observation  platform  that  could  be  guided  by  a  fiber  optic  data  link. 
Table  5. 1  presents  the  basic  aircraft  specifications. 

The  aircraft  is  somewhat  unconventional  in  that  the  engine  is  mounted  in  a  nacelle 
tractor  style  above  the  fuselage  and  wing.  The  aft  fuselage  consists  of  a  1.75  in.  diameter 
aluminum  tube  which  connects  the  tail  surfaces  to  the  main  fuselage.  Figure  5.1  displays  a 
three  view  drawing  of  the  FROG  UAV. 
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PARAMETER 

MEASUREMENT/UNITS 

Length 

8  .125  ft 

97.5  m 

Height 

1.75  ft 

21  in 

Weight 

67.7  lbs 

Power  Plant 

12Hp/  2  Cycle 

Wing  Airfoil 

NACA  2415 

Horiz.  Stab.  Airfoil 

NACA  0006  (Approx.) 

Sw(Stef) 

17.57  ft^ 

2530  in^ 

St 

3.174  ft^ 

457.1  in^ 

Sv 

0.9818  ft^ 

141.4  in^ 

C 

1.66  ft 

20  in 

Ct 

0.958  ft 

11.5  in 

K 

10.54  ft 

126.5  in 

bt 

3.313  ft 

39.75  in 

by 

1.25  ft 

15.0  in 

It 

4.44  ft 

53.25  in 

Iv  ■ 

4.44  ft 

53.25  in 

ARw 

6.32 

AR^ 

3.46 

AR, 

1.59 

Vh 

0.49 

Vv 

0.02 

Table  5.1  FROG  UAV  Characteristics,  after  Ref.  [1]. 
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Figure  5.1  FROG  UAV  Three- View  Drawing. 
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The  FROG  UAV,  as  operated  by  NFS,  is  equipped  with  airspeed,  angle-of-attack, 
altitude  and  control  surface  sensors.  In  addition,  a  nndniature  Inertial  Measurement  Unit 
(IMU)  captures  aircraft  attitude,  acceleration  and  body  rates.  Data  is  down  linked  to  a 
mobile  SGI  workstation  through  a  spread  spectrum  modem.  Onboard  GPS  provides 
differential  GPS  navigation  capability  with  the  ground  station  used  as  a  reference.  The 
aircraft  can  be  flown  by  conventional  radio  control  or  by  up-linking  flight  control 
commands  from  the  computer  workstation. 

Current  flight  control  development  revolves  around  the  cruise  trim  point  of  60 
m.p.h.  or  88  ft/s.  This  flight  condition  is  selected  for  the  development  of  stability 
derivative  data.  Table  5.2  lists  the  aircraft  parameters  for  the  trim  flight  condition. 


PARAMETER 

MEASUREMENT 

UNITS 

Weight 

67.73 

lbs 

IXX 

12.52 

lYY 

8.43 

slug-ft^ 

IZZ 

18.55 

slug-ft^ 

Airspeed 

60/88 

mph  and  ft/s 

Altitude 

800 

ftMSL 

Air  Density 

0.002327 

Slug/ft3 

Center  of  Gravity 

34.5% 

MAC 

trim 

0.4295 

n/a 

^  trim  (est) 

-1.3 

d^rees 

^  Etrim 

5.1 

d^rees 

Table  5.2  FROG  UAV  Trim  Flight  Condition,  after  Ref.[l]. 


C.  FROG  UAV  MODELING 

1.  General 

LOFTSMAN  is  utilized  for  the  creation  of  all  CMARC  input  file  patches  except 
for  wing  tips.  In  some  cases,  CMARC’s  more  efficient  built-in  capability  to  model 
standard  NACA  4-digit  wing  surfaces  could  have  been  used.  However,  future  studies  will 
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require  flight  control  surface  patches  meshed  with  LOFTSMAN.  Therefore,  with  growth 
provisions  in  mind,  all  patches  were  created  with  LOFTSMAN  from  the  start.  Figure  5.2 
displays  the  complete  FROG  UAV  model  with  all  patches  and  wakes  activated. 


Figures.!  FROG  UAV  Panel  Code  Model, 


Some  assumptions  are  made  to  simplify  the  modeling  process.  First,  the  horizontal 
and  vertical  stabilizers  are  modeled  with  a  NACA  0006  section.  The  actual  surfaces  are 
constructed  with  a  flat  section,  rounded  at  the  nose  and  tapered  starting  at  the  control 
surface  hinge  line  to  a  sharp  trailing  edge.  The  NACA  0006  provides  a  close 
approximation  and  allows  the  use  of  LOFTMAN’s  built-in  wing  lofting  capability.  For  a 
potential  flow  solution,  this  simplification  is  considered  minor. 

A  second  simplification  is  made  regarding  the  vertical  stabilizer’s  tip  rib 
orientation.  The  actual  rib  is  canted  down  5°  with  respect  to  the  longitudinal  waterline. 
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LOFTSMAN  will  only  model  a  chord  line  that  is  parallel  to  the  waterline  (constant  BL). 
The  vertical  tail  tip  rib  is  modeled  with  a  constant  BL,  but  the  span  is  adjusted  to  maintain 
the  same  overall  surface  area. 

Finally,  there  is  no  attempt  to  model  the  tricycle  landing  gear  struts  or  wheel 
assemblies.  The  landing  gear  components  do  not  contribute  significantly  to  the 
aerodynamic  stability  derivatives.  However,  they  certainly  need  to  be  taken  into  account 
when  measuring  moments  of  inertia  for  a  dynamic  model. 

2.  Modeling  Coordinate  System 

The  model  is  developed  using  a  coordinate  system  selected  to  simplify  fuselage 
measurements.  The  +x-axis  starts  even  with  the  nose  and  runs  aft  along  the  bottom  of  the 
fuselage,  parallel  with  the  tail  boom.  The  bottom  of  the  fuselage  is  used  as  the  waterline 
with  +z-axis  in  the  up  direction.  This  allows  for  easy  vertical  measurements  when  the 
aircraft  is  placed  flat  on  a  horizontal  surface.  The  +y-axis  runs  from  centerline  outboard 
parallel  to  the  right  wing.  Figure  5.3,  which  displays  static-pressure  source  and  alpha  vane 
locations,  also  shows  the  location  and  origin  of  the  modeling  coordinate  system. 

3.  LOFTSMAN  Patches 

LOFTSMAN  is  used  to  generate  all  the  model  patches  except  for  wing  tips. 
CMARC’s  built-in  capability  is  used  to  create  wing  tip  patches.  Appendix  H  contains 
listings  of  all  the  LOFTSMAN  input  files.  Once  a  surface  is  meshed,  the  mesh  is  saved  to 
a  file  as  a  CMARC/PMARC  patch.  The  resulting  text  file  is  then  opened,  and  the  text  is 
copied  and  pasted  with  any  text  editor  into  the  patch  definition  section  of  the  CMARC 
input  file.  LOFTSMAN  patch  files  are  not  listed  because  they  are  redundant  with  the 
patches  in  the  final  CMARC  input  file  listed  in  Appendix  I. 

When  saving  a  patch,  LOFTSMAN  automatically  takes  care  of  all  CMARC  input 
file  formatting  except  for  the  TNODS  patch  continuation  or  final  patch  toggle.  A  patch, 
as  formatted  by  LOFTSMAN,  assumes  additional  patches  will  follow  in  the  CMARC  input 
file.  Therefore,  the  last  segment’s  TNODS  variable  is  set  TNODS=3.  When  the  patch  is 
the  last  patch  in  the  input  file,  the  TNODS  variable  must  be  manually  set  to  TNODS=5.  If 
CMARC  hangs  up  while  reading  in  geometry  information,  most  likely  TNODS=5  is 
missing  on  the  last  patch. 
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a.  Fuselage  Model 

The  fuselage  is  lofted  as  a  B-type  body.  A  B-type  body  is  used  when  major 
portions  of  the  fuselage  have  a  circular  or  oval  cross  section.  The  input  file  is  listed  in 
Appendix  H.  Only  the  right  side  is  meshed,  with  a  symmetric  left  side  created  by  toggling 
the  DPATSYM  variable  to  IPATSYM=1.  LOFTSMAN  assumes  that  B-type  bodies 
converge  to  a  specific  point  at  the  fore  and  aft  ends.  The  flat  aft  fuselage  face  does  not 
provide  this  single  point.  A  slight  modification  was  made  to  the  aft  face  to  allow 
automatic  meshing  as  a  B-type  body.  The  center  of  the  aft  face  is  extended  very  slightly, 
approximately  1/8  inch,  to  provide  a  convergence  point  for  the  final  rear  triangular  panels. 
This  small  deviation  is  assumed  not  affect  the  aerodynamic  fidelity  of  the  model  for  a 
potential  flow  solution. 

The  right  side  was  originally  meshed  separately  from  the  wing  as  a  20  x  20 
panel  patch.  This  created  a  low  order  fit  when  the  wing  patch  was  butted  to  the  side  of 
the  fuselage,  resulting  in  overlapping  panels.  A  final  mesh  was  created  that  flowed  around 
the  wing  root  and  fuselage  intersection  for  a  high  order  fit.  All  the  fuselage  panels  at  the 
wing  root  join  with  the  adjacent  wing  panels.  This  mesh  requires  that  the  fuselage  be 
broken  up  into  six  separate  panels  per  side.  They  are  the  nose  patch,  the  forward 
transition  patch,  the  top  and  bottom  wing  root  patches,  the  aft  transition  patch  and  finally 
the  rear  fuselage  patch.  Some  manual  editing  is  required  to  straighten  out  panels  on  the 
upper  fuselage  patch.  When  the  six  patches  are  added  together,  the  final  configuration  is 
modeled  with  a  44x15  panel  patch. 

b.  Main  Wing  Patch 

The  NACA  2415  wing  is  created  with  four  separate  patches  to  allow  the 
addition  of  an  aileron  mesh  at  a  later  date.  CMARC  comes  with  a  broad  selection  of 
“*.SD”  airfoil  template  files  that  are  automatically  loaded  during  installation.  The 
“NACA2415.SD”  file  is  used  for  this  model.  The  inboard  patch  runs  from  the  wing  root, 
past  the  flaps,  to  the  start  of  the  aileron.  The  mid  patch  covers  the  portion  of  the  wing 
spanned  by  the  aileron.  The  outboard  patch  creates  the  tapered  wing  extension.  Finally,  a 
semi-circular  wing  tip  patch  is  added  in  the  input  file  using  CMARC’s  built-in  wing  tip 
functionality.  The  wing  is  set  to  a  4.5°  incidence  in  the  LOFTSMAN  input  file. 
Alternatively,  the  patch  could  be  created  with  zero  incidence  and  then  the  patch 
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coordinate  system  could  be  rotated  in  the  CMARC  input  file.  Together,  the  four  wing 
patches  add  to  make  a  20  x  30  panel  wing  model. 


c.  Horizontal  Stabilizer  Patch 

The  horizontal  stabilizer  patch  is  created  with  a  single  10  x  22  mesh  using 
the  “NACA0006.SD”  airfoil  template.  No  special  modifications  are  required.  A  tip  patch 
is  not  added  because  some  of  the  resulting  panels  would  be  too  small.  In  particular,  the 
triangular  panels  closing  out  the  aft  end  of  the  tip  are  too  small  in  proportion  to  the  other 
panels.  An  attempt  was  made  to  model  horizontal  and  vertical  stabilizer  wing  tips,  but  the 
model  will  not  converge  with  them.  Leaving  off  tip  patches  will  not  significantly  influence 
results  according  to  the  CMARC  User's  Guide  [Ref  2], 

d.  Vertical  Stabilizer  Patch 

The  vertical  stabilizer  patch  is  created  with  a  single  8x18  mesh  using  the 
“NACA0006.SD”  airfoil  template.  The  LOFTSMAN  input  file  is  different  in  that  a 
vertical  wing  surface  requires  a  modification  to  the  rib  axis.  The  rib  axis  must  be  specified 
with  an  x-axis  rotation  of  90®,  a  y-axis  rotation  of  0®  and  an  unspecified  (999.0)  z-axis 
rotation.  No  symmetry  is  selected  for  the  vertical  stabilizer  because  the  patch  is  already 
symmetric  about  the  y=0  plane.  As  with  the  horizontal  stabilizer,  a  tip  patch  is  not  added 
because  some  of  the  resulting  panels  would  be  too  small. 

e  Tail  Boom  Patch 

The  tail  boom  patch  is  created  as  a  single  12  x  10  mesh  using  a  B-type 
body.  Again,  only  the  right  side  is  meshed  due  to  symmetry.  The  LOFTSMAN  input  file 
requires  modifications  at  both  ends  in  a  similar  fashion  to  the  aft  fuselage.  A  single  point 
is  added  to  allow  convergence  of  the  triangular  panels  at  either  end.  With  this  point,  the 
tail  boom  has  the  appearance  of  being  tapered  at  both  ends.  The  point  is  then  manually 
edited  out  in  the  CMARC  input  file  by  replacing  the  “x”  coordinate  of  the  beginning  and 
ending  section  panels  with  the  correct  value.  In  most  cases,  the  tail  boom  is  left  out  of 
solution  to  aid  in  convergence.  This  is  due  to  the  small  overlapping  panels  at  the  fuselage 
tail  boom  junction.  Being  a  slender,  round  tube  directly  in  the  fuselage  slip  stream,  the  tail 
boom  should  have  little  influence  on  the  stability  derivatives. 
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f.  Engine  Pod  Patch 

The  engine  pod  patch,  or  nacelle,  is  created  as  a  single  15  x  10  mesh  using 
a  B-type  body.  Only  the  right  side  is  meshed  due  to  symmetry.  The  prop  spinner  is  an 
integral  part  of  the  patch.  No  attempt  is  made  to  model  the  prop,  engine  heads  or  exhaust 
system. 

g.  Engine  Pylon  Patch 

The  engine  pylon  patch  is  modeled  with  a  single  15x10  mesh  using  an  A- 
type  body.  A-type  bodies  are  used  to  model  surfaces  similar  to  boat  hulls  with  cornered 
surfaces  or  sharp  chines.  In  addition,  A-type  bodies  do  not  require  the  body  to  be 
completely  enclosed.  As  a  result,  an  A-body  was  selected  to  model  just  the  sides  of  the 
pylon.  Only  the  right  side  is  meshed  due  to  symmetry.  A  low  order  fit  is  achieved  with 
the  adjacent  fuselage  and  engine  pod  panels.  This  results  in  questionable  pressure 
distributions.  As  a  result,  the  pylon  patch  was  turned  off  for  most  configurations.  A 
future  attempt  "will  be  made  to  create  a  high  order  fit  between  the  other  patches.  This  will 
probably  require  manual  editing  of  the  intersecting  patches. 

4.  Common  CMARC  Input  File  Errors 

The  patches  created  in  LOFTSMAN  are  assembled  into  a  single  CMARC  input  file 
with  any  text  editor.  A  default  minimum  input  file  comes  with  CMARC  or  any  old  file 
may  be  modified.  There  are  many  errors  that  will  cause  CMARC  to  hang  up  without  an 
error  message.  The  two  most  common  errors  are  forgetting  to  designate  the  last  patch 

and  incorrectly  numbering  the  wake  patches. 

The  last  patch  must  be  designated  by  including  a  TNODS=5  setting  in  the  last 
section  of  the  last  patch.  If  it  is  not  included,  CMARC  hangs  up  when  reading  in  the 
geometry.  In  a  similar  manner,  the  last  wake  must  be  designated  with  a  NODEW=5 
setting.  If  the  last  wake  is  not  designated,  CMARC  hangs  up  while  reading  in  the  wake 
information. 

Another  common  error  involves  incorrect  wake  to  patch  number  association. 
Patch  numbering  changes  whenever  patches  are  disabled  or  reordered.  The  KWPACH 
field  for  each  wake  definition  must  be  checked  to  make  sure  it  reflects  the  current  patch 
numbering. 
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D.  STATIC-PRESSURE  SOURCE  AND  YAW  VANE  CORRECTIONS 

THROUGH  OFF-BODY  FLOW  ANALYSIS 

CMARC  is  ideally  suited  for  off-body  flow  analysis.  Off-body  streamlines  may  be 
placed  through  a  point  anywhere  in  the  flow  field.  CMARC  will  then  follow  the 
streamline  up  and  downstream  the  distance  designated  in  the  input  file.  This  is  particularly 
useful  for  flow  visualization.  In  addition,  CMARC  calculates  pressure  coefficient  and 
velocity  at  each  point  along  the  streamline.  For  this  study,  two  streamlines  are  placed 
through  the  locations  of  the  static-pressure  source  and  alpha  probe  locations.  Pressure 
coefficient  is  used  to  quantify  static  source  position  error  and  velocity  is  used  to  calculate 
alpha  probe  position  error  as  a  function  of  FROG  UAV  angle-of-attack.  Both  static 
pressure  and  AO  A  are  digitized  for  down  link  to  the  ground  station  allowing  the  values  to 
be  easily  corrected.  Either  a  look-up  table  or  curve  fit  correction  can  be  applied 
subsequent  to  being  passed  to  the  flight  control  routines. 

1.  Description  of  the  FROG  UAV  Pitot-Static  and  AOA  Systems 

The  pitot-static  system  and  angle-of-attack  probe  share  a  common  flight  test  boom 
extending  from  the  nose  of  the  UAV.  The  boom  contains  both  the  total  and  static 
pressure  ports.  Figure  5.3  depicts  the  general  dimensions  of  the  flight  test  boom 
installation  and  the  modeling  coordinate  system. 

2.  Modeling  Off-Body  Streamlines 

Streamlines  are  placed  at  the  two  locations  indicated  in  Figure  5.3  which 
correspond  to  the  static  source  and  alpha  probe  locations.  Two  off  body  streamlines  were 
activated  in  CMARC  by  setting  NSTLIN=2  in  the  &SLIN1  line.  Only  a  short  distance  of 
2  inches  is  selected  up  and  downstream  in  the  SU  and  SD  fields  to  reduce  the  size  of  the 
output  file.  Figure  5.4  is  a  POSTMARC  rendering  of  the  two  off-body  streamlines  used 
for  sensor  corrections.  With  the  model  at  at=0°,  notice  that  the  streamline  is  curving  up  at 
the  angle-of-attack  vane  location  6.5  inches  in  front  of  the  aircraft  nose. 
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Figure  5.3  Diagram  of  the  FROG  UAV  Pitot-Static  and  AOA  Systems. 
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Figure  S.4  FROG  UAV  Off-Body  Streamline  visualization  with  POSTMARC 

(at=10°). 

3.  Analysis  of  Static  Source  Position  Errors 

In  general,  the  position  error  pressure  coefficient,  ACp  pc  or  APp/qc,  is  a  function  of 
ffeestream  Mach  number  and  angle-of-attack  provided  that  the  static  source  is  located 
outside  of  a  thick  boundary  layer  and  sideslip  is  minimized  [Ref  13],  In  the  case  of  the 
FROG  UAV  with  incompressible  flow,  APp/qc  becomes  a  function  of  angle-of-attack  only. 
As  a  result,  the  corrections  can  be  simply  defined  as  a  function  of  measured  angle-of- 
attack. 

A  DOS  batch  file  was  executed  to  step  the  CMARC  model  through  angles-of- 
attack  ranging  from  -8’’  to  20°.  The  batch  file  incremented  the  angle-of-attack  using 
CMARC’s  command  line  override  feature.  In  addition,  a  new  output  file  name  was 
designated  for  each  angle-of-attack.  Position  error  pressure  coefficient  is  then  read  from 
the  off-body  streamline  listing  of  the  output  file  at  the  location  corresponding  to  the  static 


source.  Table  5.3  lists  the  values  of  APp/qc  calculated  from  CMARC  data.  Figure  5.5 
displays  ACp  pc  as  a  function  of  indicated  angle-of-attack.  The  second  order  influence  of 
angle-of-attack  is  clear  with  the  second  order  curve  fitting  tightly  through  the  data  points. 
Of  note,  the  error  is  relatively  constant  for  a  ±8°  band  around  trim  angle-of-attack.  For 
incompressible  flow,  position  error  pressure  coefficient  is  independent  of  airspeed  and 
altitude. 

Position  error  pressure  coefficient  can  be  turned  into  position  corrections  for 
airspeed  and  altitude.  The  following  relations  were  developed  which  assume  small  errors 
and  incompressible  flow: 

VjACp 

AVpc= — and  AVpc=Vc-Vi  5.1 

AVncVi 

AHpc=— ^  and  AHpc=Hc-Hi  5.2 

^  c^stdgo 

Where: 

AHpc  is  the  altitude  position  correction. 

AVpc  is  the  velocity  position  correction. 

APn  rc  ■ 

ACn  =  — ^  or  position  error  pressure  coefficient. 

^  qc 

Ostd  is  standard  day  density  ratio, 
go  is  the  gravitational  constant. 

Table  5.3  displays  corrections  calculated  for  both  airspeed  and  altitude  at  the 
FROG  UAY  trim  condition  of  88  ft/s  and  800  ft  MSL.  The  corrections  are  added  to  the 
indicated  value  to  obtain  the  corrected  value.  Figures  5.6  and  5.7  display  the  corrections 
as  a  function  of  indicated  angle-of-attack.  Again,  a  second  order  curve  fits  nicely  through 
the  data  points.  Equations  5.1  and  5.2  can  be  used  to  implement  a  correction  algorithm 
based  on  airspeed  and  altitude. 
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Table  5.3  Position  Error  Corrections  for  the  NPS  FROG  UAV  at  V=88  ft/s  and 
H=800  ft  MSL.  Derived  from  CMARC  Panel  Code  Off-Body  Flow 
Field  Analysis. 


Source:  CMARC  Panel  Code 
Trim  Condition;  V=88  ft/s  H=800  ft  MSL 


Figure  5.5  Position  Error  Pressure  Coefficient,  ACp  pc,  for  the  NPS  FROG  UAV. 

Derived  from  CMARC  Panel  Code  Off-Body  Flow  Field  Analysis. 


I 
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AHpc  (ft) 


Source:  CMARC  Panel  Code 
T rim  Condition;  V=88  ft/s  H=800  ft  MSL 
Flight  Test  Boom  Static  Source  (8.25"  forward  of  nose) 


AOA  Indicated  (deg) 


Figure  5.6  Altitude  Position  Error,  AHpc,  for  the  NPS  FROG  UAV  at  V=88  ft/s 
and  H=800  ft  MSL.  Derived  from  CMARC  Panel  Code  Off-Body 
Flow  Field  Analysis. 
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Source:  CMARC  Panel  Code 
Trim  Condition;  V=88  ft/s  H=800  ft  MSL 
Flight  Test  Boom  Static  Source  (8.25"  forward  of  nose) 


AOA  Indicated  (deg) 


Figure  5.7  Airspeed  Position  Error,  AVpc,  for  the  NFS  FROG  UAV  at  V=88  ft/s 
and  H=800  ft  MSL.  Derived  from  CMARC  Panel  Code  Off-Body 
Flow  Field  Analysis. 
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4.  Analysis  of  Alpha  Vane  Position  Error 


Local  flow  field  velocity  is  extracted  fi-om  the  ofF-body  streamline  listing  to  obtain 
local  angle-of-attack.  The  alpha  vane  is  assumed  to  capture  the  x-z  component  of  the 
local  velocity  field  and  ignore  cross  flow  in  the  y  direction.  Flow  field  velocity  is  turned 
into  indicated  angle-of-attack  and  angle-of-attack  position  correction  with  the  following 
equations: 


ai°=  atan 


^  z 


Wxy 


*180  ^ 

* -  degrees 

% 


5.3 


Attpc  =at  -aj  degrees  5.4 

A  DOS  batch  file  is  executed  to  step  the  CMARC  model,  with  an  ofF-body 
streamline  located  at  the  vane  position,  through  angles-of-attack  ranging  fi-om  -8°  to  20°. 
Local  velocity  components  are  then  read  from  the  location  corresponding  to  the  alpha 
vane.  Table  5.4  lists  the  values  of  Aopc  calculated  from  CMARC  data.  Figure  5.8  displays 
Aa  pc  as  a  function  of  indicated  angle-of-attack.  Linear  and  second  order  curve  fit 
equations  are  also  indicated  on  Figure  5.8.  Angle-of-attack  correction  is  fairly  linear 
through  the  FROG  operating  envelope,  with  approximately  -1.25  degrees  of  position  error 
at  the  FROG  cruise  trim  condition.  The  corrections  apply  at  all  incompressible  airspeeds 
and  all  altitudes. 

5.  Summary  of  Off-Body  Flow  Field  Analysis 

CMARC  proved  useful  for  both  static-pressure  source  and  alpha  vane  position 
corrections.  Measured  data  may  be  corrected  using  look-up  tables  with  the  values  in 
Table  5.3  and  5.4  or  by  using  the  curve  fits  in  Figures  5.5  through  5.8.  Flight  testing  is 
recommended  for  validation  of  sensor  corrections  obtained  from  this  CMARC  ofF-body 
flow  field  analysis. 
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UAV  AOA  Velocity  at  Alpha  Vane 

aT(d^)  Vx(ft/S) 


81.60 


81.67 


.71 
.73 
.73 
81.70 


81.66 

81.58 

81.48 

81.37 

81.07 

80.67 

80.17 
79.61 
78.93 

78.18 
77.34 


Vy  (ft/s) 

1.66 

1 

.65 

1 

.64 

1 

.63 

1 

.63 

1 

.62 

1 

.61 

1 

.60 

1 

.59 

1 

.58 

1 

.57 

1 

.56 

1 

.53 

1 

.51 

1 

.48 

1 

.46 

-12.23 


-1.47 


0.28 

2.13 


3.93 

5.72 

7.51 


9.30 

11.08 

12.88 

16.43 

19.98 

23.50 

26.99. 


37.31 


AOA  Correction 
Aa=aT- tti  (d^) 

AOA  Indicated 
ai(deg) 

0.60 

-8.60 

0.08 

-6.08 

-0.75 

-2.25 

-0.97 

-1.03 

-1.20 

0.20 

-1.49 

1.49 

-1.75 

2.75 

-2.00 

4.00 

-2.25 

5.25 

-2.50 

6.50 

-2.75 

7.75 

-2.99 

8.99 

-3.46 

11.46 

-3.91 

13.91 

-4.34 

16.34 

-4.73 

18.73 

-5.11 

21.11 

-5.44 

23.44 

-5.75 

25.75 

Table  5.4  Angle-of  Attack  Vane  Position  Error  Corrections  for  the  NPS  FROG 
UAV.  Derived  from  CMARC  Panel  Code  Off-Body  Flow  Field 
Analysis. 
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Source:  CMARC  Panel  Code 
Trim  Condition:  V=88  ft/s  H=800  ft  MSL 
Flight  Test  Boom  Alpha  Vane  (6.5"  forward  of  nose) 


-12  -8  -4  0  4  8  12  16  20  24  28 

AOA  Indicated  (deg) 


Figure  5.8  Angle-of-Attack  Vane  Position  Error,  Attpc,  for  the  NFS  FROG  UAV. 

Derived  from  CMARC  Panel  Code  Off-Body  Flow  Field  Analysis. 
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E.  DEVELOPMENT  OF  BASIC  STABILITY  DERIVATIVES 


In  this  section,  CMARC  is  used  to  develop  some  of  the  basic  longitudinal  and 
lateral-directional  stability  derivatives  for  the  FROG  UAV.  The  development  effort 
focuses  on  finding  the  Clo  and  Cma  longitudinal  stability  derivatives  followed  by  the  Cyp, 
Cip  and  Cnp  lateral-directional  stability  derivatives.  Control  power  and  rate  damping 
derivatives  will  be  the  focus  of  ongoing  research. 

CMARC  contains  built-in  functionality  to  integrate  forces  and  moments  in  all  axes 
over  the  surface  of  a  body.  Forces  and  moments  are  automatically  normalized  into  non- 
dimensional  coefficients  based  on  the  mean  aerod3mamic  chord,  reference  wing  area,  semi¬ 
span  and  center  of  gravity  location  in  the  CMARC  BINP9  input  line.  Coefficients  are 
presented  in  both  ^vind  and  body  axes.  The  CMARC  model  is  run  at  two  different  angles- 
of-attack  and  one  sideslip  angle.  The  slope  of  the  force  and  moment  coefficients  is  then 
taken  to  produce  the  Ci^x  and  Cma  longitudinal  derivatives  and  the  Cyp,  Qp  and  Cnp  lateral- 
directional  derivatives. 

The  CMARC  model  must  be  analyzed  in  the  linear  slope  regions  of  a  and  P  for 
valid  results.  A  potential  flow  solution  will  not  produce  satisfactory  results  for  bodies 
with  significant  areas  of  flow  separation. 

1.  Longitudinal  Stability  Derivatives 

a.  Longitudinal  Stability  Derivative  Methods 

Three  basic  longitudinal  stability  derivatives  can  be  measured  with  just  two 
runs  of  the  CMARC  model.  The  model  is  first  analyzed  at  an  angle-of-attack 
corresponding  to  the  estimated  trim  condition.  In  this  case,  at=0°  is  selected  for  the  first 
run.  A  second  CMARC  run  is  conducted  with  angle-of  attack  incremented  one  or  two 
degrees.  Cl  and  Cm  are  then  extracted  manually  from  the  data  files.  The  slope  of  Cl  and 
Cm  versus  angle-of-attack  provide  the  CLa  and  Cma  longitudinal  derivatives.  For  this 
study,  several  angles-of-attack  were  analyzed  to  check  consistency  of  the  slope.  In 
addition,  auim  is  calculated  fi'om  the  lift  curve  slope  and  trim  lift  coefficient.  Equations  5.5 
through  5.7  are  used  for  these  calculations.  For  the  longitudinal  analysis,  only  half  the 
model  is  analyzed.  The  symmetric  calculation  mode  is  selected  by  setting  both 
RSYM=0.0  and  IPATSYM=0  in  the  CMARC  input  file. 


75 


per  radian 


5.5 


C  _fcL2-CLi)^180 
(a2-ai) 


(^m2  ^mi)  180 
C  ma  =  “7 - P  ~  per  radian 


trim 


:a°i  + 


*180 
* _ 

n 


degrees 


5.6 


5.7 


Several  FROG  UAV  model  configurations  were  analyzed  in  a  build-up 
approach  to  check  results  against  classical  calculations  and  flight  test  data.  Figure  5.9 
shows  the  simplified  CMARC  models.  First,  just  the  wing  and  horizontal  tail  were 
considered.  The  patches  for  all  other  surfaces  and  wakes  were  turned  off  and  the  wing 
root  was  extended  to  centerline.  The  FROG  fuselage  was  then  analyzed  separately  and 
the  results  were  added  to  the  simplified  wing  and  horizontal  tail  combination.  Next,  the 
original  butted  (low  order  fit)  wing/fuselage  and  horizontal  tail  were  considered.  Finally, 
the  blended  wing/fuselage  and  horizontal  tail  were  analyzed.  Values  of  Cta  and  for 
these  four  configurations  are  presented  in  Table  5.5. 

Classical  design  calculations  are  also  performed  to  estimate  Cma  for 
comparison  to  CMARC  results.  Equation  5.8  is  used  for  the  calculation  of  Cma-  In 
classical  design,  the  horizontal  tail  downwash  derivative,  de/da,  is  generally  selected  from 
empirical  data.  Using  a  taper  ratio  of  TR=1;1  and  aspect  ratio  of  AR=6,  d8/da=0.4  is 
selected  from  empirical  charts  in  Ref  [14]  for  the  FROG  UAV  configuration.  A  few  other 
values  of  the  horizontal  tail  downwash  derivative,  de/da,  are  selected  to  see  how  well 
CMARC  models  downwash  effects.  Classical  design  estimates  of  Cma  for  values  of  de/da 
ranging  from  0  to  0.4  are  presented  in  Table  5.5  for  comparison  with  CMARC  results. 
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Figure  5.9  Simplified  CMARC  Models  of  the  FROG  UAV  for  the  Determination 
of  Longitudinal  Stability  Derivatives. 
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Flight  test  data  for  the  short  period  and  phugoid  modes  were  used  for 
longitudinal  parameter  estimation.  Values  for  Cux  and  based  on  preliminary 
parameter  estimation  work  by  Engdahl  [pending  publication]  are  presented  in  Table  5.5. 
Caution  is  advised  against  making  definitive  comparisons  until  the  work  is  published. 


METHOD 

CONFIGURATION* 

LONGITU 

Ottrim 

(aeg) 

DINAL  PAR 

^La 

(per  raa; 

i\METERS 

Cma 
(per  raa; 

CMARC 

Panel  Code 

Wing/Horiz  Tail 

-0.87 

4.86 

-0.835 

Wing/Horiz  Tail  +  Fuselage 

4.78 

-0.608 

Blended  Wing-Fuselage/Horiz  Tail 

-0.01 

4.72 

-1.105 

Butted  Wing-Fuselage/Horiz  Tail 

-0.8 

5.37 

-1.348 

Classical 

Design^ 

Wing/Horiz  Tail  -  5s/5a=0 

-0.78 

4.89 

-1.50 

Wing/Horiz  Tail  -  6s/5a=0.25 

-0.81 

4.85 

-1.00 

Wing/Horiz  Tail  -  58/5a=0.35 

-0.82 

4.83 

-0.80 

Wing/Horiz  Tail  -  6E/8a=0.40 

-0.82 

4.82 

-0.70 

Parameter 

Estimation"* 

Flying  Aircraft 

n/a 

4.09 

-0.42 

NOTES:  1)  CGx=34.5%  M.A.C.  /  062=8.6"  from  bottom  of  fuselage. 

2)  Zero  lift  wing  incidence  is  +6.5'’  from  the  longitudinal  reference  line. 

3)  Classical  Design  after  Ref.  [14]. 

4)  Unpublished  parameter  estimation  from  flight  test  data  by  Engdahl. 


Table  5.5  Comparison  of  FROG  UAV  Longitudinal  Stability  Derivatives. 

b.  Analysis  of  Longitudinal  Stability  Data 

The  first  three  configurations  in  Table  5.5  produce  good  results  for  Ctaand 
reasonable  values  for  Cma-  However,  the  fourth  configuration,  the  butted  wing  root  and 
fuselage,  produces  excessively  large  values  of  both  Clo  and  Cmo-  This  configuration 
should  be  avoided  in  future  models.  It  is  recommended  that  CMARC  model  developers 
spend  the  time  up  front  to  produce  the  higher  fidelity  model  from  the  start. 

The  values  produced  for  Cma  from  CMARC  are  somewhat  high  when 
compared  with  to  the  classical  design  calculations.  Clearly,  some  downwash  is  sensed  by 
the  horizontal  tail  in  the  CMARC  analysis  because  all  values  for  Cma  are  considerably  less 
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than  the  classical  calculation  with  d8/da=0.  Still,  high  values  compared  to  flight  test  data 
indicates  that  CMARC  has  a  difficult  time  capturing  the  complete  ds/da  downwash  effect. 
This  could  be  due  to  the  requirement  to  select  rigid  wakes  to  prevent  the  wing  wake  from 
impacting  the  horizontal  tail,  A  more  careful  wake  definition  may  help  capture  the  tail 
downwash  derivative  with  more  fidelity.  A  study  by  Walden  et  al.  [Ref.  15]  studied  wake 
turbulence  by  modeling  an  aircraft  flying  in  trail  of  a  wake  generating  wing.  A  horizontal 
tail  trailing  the  main  wing  is  a  similar  configuration.  The  study  found  that  a  streamline- 
based  wake  is  the  best  method  for  modeling  downwash  effects.  This  wake  definition 
should  be  investigated  for  modeling  the  Cma  derivative.  Of  note,  the  wake  diffusion 
process  is  neglected  ia  a  potential  flow  analysis. 

In  summary,  CMARC  produced  accurate  values  for  Cta  and  slightly  high 
values  of  Cma-  Difficulties  were  encountered  trying  to  model  the  horizontal  tail  downwash 
derivative.  A  more  careful  study  of  the  effects  of  wing  wake  placement  on  the  downwash 
derivative  is  recommended. 

2.  Lateral  Directional  Stability  Derivatives 

a.  Lateral-Directional  Stability  Derivative  Methods 

Development  of  the  lateral-directional  stability  derivatives  is  more  straight 
forward  than  for  the  longitudinal  derivatives  because  the  vertical  tail  sidewash  angle  plays 
a  lesser  role.  However,  both  sides  must  be  modeled  by  setting  both  RSYM=1.0  and 
IPATSYM=1.  This  creates  symmetric  patches  around  the  y=0  plane  allowing  CMARC  to 
perform  asymmetric  calculations  around  the  entire  body  and  significantly  increases 
processing  times. 

For  the  lateral-directional  axis,  the  aircraft  is  modeled  with  the  blended 
wing  and  fuselage  in  combination  with  the  vertical  and  horizontal  stabilizers  as  shown  in 
Figure  5.10.  The  engine  nacelle  and  pylon  are  left  off  because  their  wakes  impact  the 
vertical  tail.  In  addition,  the  pylon/fuselage  and  pylon/nacelle  junctions  were  meshed  with 
a  low  order,  butted  fit.  This  type  of  junction  was  found  to  produce  poor  results  during  the 
longitudinal  stability  study. 

The  model  is  first  checked  for  lateral  directional  balance  at  zero  yaw  angle. 
The  side  force,  rolling  and  yawing  coefficients  should  be  zero  when  a  trial  run  is 
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performed  at  zero  yaw  angle.  If  lateral-directional  forces  or  moments  are  present,  the 
model  and  wake  geometry  should  be  checked  for  symmetry. 


Figure  5.10  Simplified  CMARC  Model  of  the  FROG  UAV  for  the  Determination 
of  Lateral-Directional  Stability  Derivatives. 

Next,  a  single  CMARC  run  is  performed  with  a=  atnm  and  yaw  angle  set  to 
one  degree.  The  lateral-directional  derivatives,  Cyp,  Cip  and  Cnp,  are  then  obtained  directly 
with  equations  5.9  through  5.11; 


Cy  *  180 

Cvo  = - * - per  radian 

Ap°  %  ^ 

5.9 

Cl  *180 

Ci„  - - *  per  radian 

5.10 

Cn*180 

C„„  -  *  per  radian 

"P  Ap°  %  ^ 

5.11 

It  should  be  noted  that  the  stability  axis  as  modeled  (x-aft  and  z-up)  differs 
from  the  standard  flight  dynamics  stability  axis.  Care  must  be  taken  to  reverse  the  signs  of 
the  appropriate  coefficients  to  convert  from  a  CMARC  model’s  stability  axis  into  the  flight 
dynamics  stability  axes 

b.  Analysis  of  Lateral-Directional  Stability  Data 

Lateral-directional  stability  derivatives  obtained  from  CMARC  are 
presented  in  Table  5.6.  For  comparison  three  other  sets  of  data  are  also  presented.  The 
first  comes  from  the  classical  analysis  presented  by  Papageorgiou  in  Ref  [1].  The  second 
set  comes  from  estimates  based  on  data  recorded  from  flight  test  static  sideslip  maneuvers, 
also  published  in  Ref  [1].  The  third  set  comes  from  parameter  estimation  by  Engdahl 
based  on  dynamic  flight  test  data.  It  is  unpublished  and  should  be  considered  preliminary. 
In  all  cases,  the  CMARC  lateral-directional  stability  derivatives  produce  a  closer  match  to 
flight  test  data  than  those  derived  from  classical  methods.  It  is  concluded  that  CMARC  is 
a  good  tool  for  lateral-directional  stability  analysis. 

3.  Summary  of  CMARC  Stability  Derivative  Analysis 

In  summary,  for  the  longitudinal  axis,  CMARC  produces  accurate  values  for  (Xtnm 
and  CLa  and  slightly  high  values  of  C„a.  Difficulties  may  be  encountered  while  trying  to 
model  the  horizontal  tail  downwash  derivative.  A  more  careful  study  of  the  effects  wing 
wake  placement  on  the  downwash  derivative  is  recommended.  Specifically,  modeling 
should  include  streamline-based  wake  placement  techniques  [Ref  15].  Analysis  of  the 
lateral-directional  axis  proves  more  straightforward.  Lateral-directional  derivatives  from 
CMARC  for  Cyp,  Qp  and  Cnp  provide  a  closer  match  to  flight  test  data  than  the  classical 
estimates.  The  engine  nacelle  and  pylon  should  be  re-meshed  and  included  in  future 
studies. 
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Overall,  CMARC  derived  stability  derivatives  are  sufficiently  accurate  for 
entry  into  an  initial  aerodynamic  model.  Adjustments  through  analysis  of  flight  test  data 
will  still  be  required.  Future  CMARC  studies  should  concentrate  on  developing  the  rate 
damping  and  control  power  derivatives. 


METHOD 

configuration’ 

LAT-D 

RPARAM 

ETERS 

Cyp 

(per  raaj 

C|p 

iper  rao) 

Cnp 
(per  rao) 

CMARC 

Panel  Code 

Blended  Wing-Fuselage/HorzA/ert  Tails 

-0.573 

-0.063 

0.120 

Classical 

Design^ 

Wing/FuselageA/ert  Tail 

-0.310 

-0.051 

0.058 

Flight  Test^ 

Flying  Aircraft 

-0.700 

-0.053 

0.057 

Parameter 

Estimation^ 

Flying  Aircraft 

-0.987 

-0.094 

0.176 

NOTES:  1)  CGx=34.5%  M.A.C.  /  003=8.6"  from  bottom  of  fuselage. 

2)  Classical  Design  calculations  by  Papageorgio,  from  Ref.  [1]. 

3)  Flight  test  results  from  Steady  Heading  Sideslip,  from  Ref.  [1] 

4)  Unpublished  parameter  estimation  from  flight  test  data  by  Engdahl. 


Table  5.6  Comparison  of  FROG  UAV  Lateral-Directional  Stability  Derivatives, 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


CMARC  is  a  DOS  personal  computer  hosted  panel  code  adopted  from  the  NASA 
Ames  PMARC  code.  AeroLogic,  Inc.,  created  CMARC  by  converting  PMARC 
FORTRAN  77  source  code  into  the  C  language.  Significant  memory  management  and 
command  line  enhancements  were  also  added.  CMARC  solves  for  inviscid, 
incompressible  flow  over  complex  three-dimensional  bodies.  Emphasis  in  this  study  is  first 
placed  on  verifying  CMARC  against  the  PMARC  and  Naval  Postgraduate  School 
Unsteady  Potential  Flow  (UPOT)  panel  codes.  CMARC  pressure  distributions  and 
boundary  layer  calculations  are  then  compared  to  experimental  data  for  an  inclined  prolate 
spheroid.  Finally  a  complex  three-dimensional  panel  model  of  the  Naval  Postgraduate 
School  FROG  UAV  is  developed  which  successfully  generates  static-pressure  source 
position  corrections,  alpha  vane  correction  curves  and  basic  stability  derivatives. 

CMARC  results  are  found  to  be  equivalent  to  the  NASA-Ames  PMARC  panel 
code.  As  expected,  pressure  distribution  and  boundary  layer  calculations  from  CMARC 
match  exactly  those  obtained  with  PMARC.  The  following  enhancements  are  noteworthy. 
CMARC,  hosted  on  a  Pentium  class  PC,  processes  input  files  significantly  faster  than 
PMARC  hosted  on  a  networked  SGI  Indigo^  UNIX  workstation.  CMARC’s  extensive 
command  line  frinctionality  greatly  enhances  batch  file  processing  capabilities.  On  the 
other  side,  CMARC’s  poor  error  flagging  capability  leaves  the  user  frequently  spending 
much  time  searching  for  input  file  mistakes.  Improved  input  file  error  checking  should  be 
incorporated  into  CMARC  functionality. 

CMARC  integral  boundary  layer  calculations  are  compared  to  the  two-dimensional 
finite  difference  methods  implemented  in  the  UPOT  code.  In  general,  CMARC  provides 
correct  trends  for  both  the  transition  and  separation  points.  However,  in  all  cases, 
CMARC  predicts  early  transition  and  late  flow  separation.  As  expected,  the  differences 
are  greatest  at  lower  Reynolds  numbers  where  boundary  layer  thickness  is  larger.  An 
adjustment  of  the  empirical  transition  and  separation  models  contained  in  CMARC  may 
prove  useful. 

CMARC  calculations  are  also  compared  to  wind  tunnel  data  for  a  6:1  inclined 
prolate  spheroid  model  at  10  degrees  angle-of-attack.  With  proper  wake  placement, 
CMARC  can  produce  accurate  normal  force  and  pitching  moment  coefficients.  Over  the 
three  dimensional  body,  CMARC  boundary  layer  calculations  also  predict  early  transition 
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and  late  flow  separation.  Despite  inaccuracies,  CMARC  boundary  layer  calculations 
remain  useful  when  used  as  a  design  tool  for  visualizing  the  trend  in  transition  and 
separation  points  with  configuration  changes. 

CMARC  integrated  skin  friction  forces  are  compared  to  prolate  spheroid  wind 
tunnel  data.  Normal,  axial,  and  pitching  moment  coefficients  for  skin  friction  forces  are 
underpredicted  by  CMARC,  but  remain  within  40%  of  integrated  experimental  data. 

The  LOFTSMAN  and  POSTMARC  portions  of  the  Personal  Simulation  Works 
software  suite  are  used  exclusively  for  the  pre-process  modeling  and  post-process 
visualization  of  CMARC  files.  The  LOFTSMAN  capability  to  automatically  format  and 
generate  CMARC  input  patches  is  an  enhancing  characteristic.  Functionality  should  be 
added  to  allow  the  modeling  of  wing  tip  ribs  that  are  not  parallel  to  the  aircraft  butt  line. 

POSTMARC  is  an  excellent  tool  for  visualizing  CMARC  output  files.  The 
capability  to  create  streamlines  and  perform  boundary  layer  calculations  external  to 
CMARC  is  extremely  useful.  However,  much  time  could  be  saved  if  POSTMARC 
maintained  previous  settings  and  selections  following  translations,  rotations  and  re-scaling. 
Additionally,  a  capability  to  overlay  multiple  data  types  is  desired. 

CMARC  off-body  flow  field  analysis  is  useful  for  both  static-pressure  source  and 
alpha  vane  position  corrections.  Measured  data  may  be  corrected  using  look-up  tables  or 
through  curve  fits  of  CMARC  derived  data.  Flight  testing  is  recommended  for  validation 
of  sensor  corrections  obtained  from  the  CMARC  off-body  analysis. 

For  the  longitudinal  analysis,  CMARC  produces  accurate  values  for  atnm  and  Cta 
and  slightly  high  values  of  Cma  Some  difficulties  are  encountered  trying  to  model  the 
horizontal  tail  downwash  derivative.  A  more  careful  study  of  the  effects  of  wing  wake 
placement  on  the  downwash  derivative  is  recommended. 

Analysis  of  the  lateral-directional  axis  proves  more  straightforward.  Lateral- 
directional  derivatives  from  CMARC  for  Cyp,  Cip  and  Cnp  provide  a  closer  match  to  flight 
test  data  than  classical  design  calculations.  Adjustments  through  analysis  of  flight  test  data 
may  still  be  required.  The  engine  nacelle  and  pylon  should  be  re-meshed  and  included  in 
future  studies. 

Overall,  the  CMARC  panel  code  is  found  to  be  suitable  for  aerodynamic  modeling 
of  the  Naval  Postgraduate  School  FROG  UAV.  CMARC  derived  stability  derivatives  are 
sufficiently  accurate  for  incorporation  into  an  initial  aerodynamic  model.  Future  CMARC 
studies  should  concentrate  on  the  development  of  the  rate  damping  and  control  power 
derivatives. 
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APPENDIX  A. 


DEVELOPMENT  OF  THE  MOMENTUM  INTEGRAL  EQUATION 

The  CMARC  and  PM  ARC  User’s  Guides  contain  the  development  of  the 
implemented  boundary  layer  equations  starting  from  the  two-dimensional  momentum 
integral  equation.  For  completeness,  the  momentum  integral  equation  is  developed  here 
to  provide  continuity. 

The  development  of  the  momentum  boundary  layer  equations  is  outlined  by 
Young  in  Ref.  [9].  In  1904  Prandtl  first  presented  his  Boundary  Layer  Theory  based  on 
the  following  observations: 

1)  However  small  the  viscosity  of  a  fluid  ,  it  cannot  be  ignored.  At  the  surface, 
the  fluid  is  at  rest  compared  to  the  body  (no  slip  condition). 

2)  Shear  stresses  are  directly  proportional  to  the  rates  of  strain. 

3)  The  ratio  of  inertial  forces  to  viscous  forces,  or  Reynolds  number,  is  important 
in  characterizing  flow  phenomena. 

4)  The  full  non-linear  viscous  Navier-Stokes  equations  are  difficult  to  solve 
directly.  Prandtl  observed  that  simplifications  could  made  when  assuming  a 
thin  boundary  layer.  Viscosity  can  be  ignored  outside  the  boundary  layer 
allowing  the  use  of  classical  inviscid  methods. 

Thin  boundary  layer  theory  also  assumes  that  the  pressure  distribution  outside  the 
thin  boundary  layer  is  transmitted  normally  through  the  boundary  layer  to  the  surface 
without  loss.  CMARC  takes  advantage  of  this  assumption  by  neglecting  the  thickness  of 
the  boundary  layer  and  imposes  a  potential  flow  solution  over  the  surface. 

The  momentum  integral  equation  for  two-dimensional  incompressible  flow  is  the 
starting  point  for  the  boundary  layer  analysis  outlined  in  References  [2]  and  [4].  It  is 
obtained  through  the  following  total  energy  integral  analysis  as  outlined  by  Young  in 
Ref.  [9], 
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Figure  A.  1  depicts  an  incremental  portion  of  a  two-dimensional  boundary  layer. 
The  mass  flow  rate  )  across  each  side  is  given  by; 


riiAD  -  0 


dx 
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Jpudz 
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Ax  -I-  O 


(ax2) 


mBc=Pe'''h^^  and  riipc  -  riiAB  =  riiRC  continuity. 
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PeWh  = 


dx 


h 

Jpudz 

LO 


+  0(Ax) 


Figure  A.l  Elementary  boundary  layer  section  for  deriving  the  momentum 
integral  equation  for  two-dimensional  flow,  after  Ref.  [9]. 
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Similarly,  the  rate  of  momentum  transport  across  each  boundary  is  given  by; 


AD:  =0 


DC-AB:  =  — 
ax 


I  pu^dw 


BC:  =pw^-u^Ax  =  u^  — 


O(Ax') 

'h 

I  pudz 


A.2 


Ax  +  0(Ax') 


The  force  due  to  pressure  on  the  sectional  boundary  layer  element  is  given  by; 
=  -hAp  =  ^  +  o(Ax^) 


A.3 


And,  the  friction  force  exerted  by  the  wall  is; 

=  ■  Ax  A.4 

Summing  the  momentum  terms  and  equating  them  to  the  forces  while  taking  the  limit  as 
Ax  0  5delds  the  momentum  integral  equation; 


dx 


\pudz 


A.5 


It  is  more  convenient  to  express  the  relation  in  terms  of  displacement  and  momentum 
thickness  by  substituting  the  following; 

dp  du,  .  . 

-  —  = 

dx  ‘  dx 

The  momentum  integral  is  then  reduced  to; 


dx 


n 

I pu{u  - 


+  - 


dx 


J  pudz 


du^ 


d_ 

dx 


Lo 


jpu(u-u^)dz 


dx 


■  -■ r„. 


d  /  7„\  du^ 


A.  7 
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A.  8 


Substituting  H  =  S*  j 6 where  H  is  the  boundary  layer  shape  factor,  and  rearranging  after 
the  chain  rule,  the  momentum  integral  can  be  written  in  as; 


ax  u  ax  Pe  ^ 


A.9 


And  finally,  by  substituting  Cf=  — 


one  obtains  Equation  16  in  References  [2] 


and  [4]: 


dx  dx  ^  dx  2 


A.  10 


From  here,  the  CMARC  or  PMARC  guides  provide  a  detailed  development  of  the 
implemented  boundary  layer  models. 
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APPENDIX  B. 


INTEGRATION  OF  AERODYNAMIC  FORCES  OVER  THE  SURFACE  OF  A 

PROLATE  SPHEROID 

The  experimental  set-up  in  Ref.  [12]  did  not  include  measurement  of  forces. 
However,  it  was  deemed  that  the  2000+  pressure  and  500+  skin  friction  measurements 
would  be  sufficient  to  allow  the  integration  of  measurements  over  the  surface  of  the 
prolate  spheroid  for  a  good  approximation  of  total  force  and  moment  coefficients.  The 
following  technique  is  developed  to  provide  an  estimate  of  integrated  pressure  and  friction 
forces.  Symmetry  is  assumed.  Appendix  C  lists  the  entire  MATLAB  program  which 
implements  the  technique  that  follows. 

In  general  the  pressure  force  is  given  by: 

Fp  =  J|  Fnc^ ,  « is  a  unit  surface  normal  B.  1 

s 

However,  the  test  data  is  provided  discreetly  in  cylindrical  coordinates,  resulting  in  the 
following  discrete  double  summation: 

Fp  =  PnrAd^Ax  /  2a ,  where  dS=  r^^Aic  /2a  B .  2 

x/2a  <f> 

Figure  B.l  shows  a  diagram  of  the  pressure  and  skin  friction  acting  over  the 
incremental  surface  areas,  AS.  The  pressure  and  skin  fiiction  coefficients,  scalars,  are 
assumed  to  be  constant  over  the  incremental  surface. 

AF^  =  PnrA^Ax  /  2a  =  (  qCp  +  IwAS ,  B.3 

where  Cp  = - —  P  =  Cp  q  +  B.4 

q 
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Flow  Direction 


Right  Side  View 


Bottom 


Ri^t 


View  Looking  Forward 


Figure  B.l  Prolate  Spheroid  G^eometry  and  Forces 


Free  stream  pressure  (P„),  assumed  to  be  constant,  can  be  dropped  from  the 
integration  due  to  symmetry.  This  leaves  the  following  relation: 

1  360 

Fp  =  ^  ^qCpnrA^/Sx  1 2a  B.5 

xl2a=0  ^=0 

Like'wise,  skin  friction  can  also  be  integrated  using  the  following  relations; 

AFsf  =  qCj^vAS  =  qCf^vrA^Ax  /  2a ,  where  v  is  a  unit  velocity  vector  B.6 
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^SF=  Z  ^rvA^Ax  I  2a 


xf2a=0  ^=0 


The  pressure  force  coefficients,  normalized  by  S  =  7d>^  and  c-2b,  are  yielded  by 
discretely  integrating  the  following  equations  in  a  cylindrical  coordinate  system: 


Q,  =  =  2Z  Z  -  i^Sprn)  ■  kA<f>JAx,  /  2a 

q<io^  >=i  j=i 


,  4  =  2Z  Z  ■  i  Ax/  /  2a 

'  <7=0*^  <=i  y=i 


Xyf  mn  _  - 

Cm,  =  =  2Z Z feooCp/ff) \xJ2a-k-zJ2a- /)]A(zJ^. Ax,.  / 2a 

1=1  y=i 


Where  the  surface  unit  normal  is  given  by: 


Unit  Normal:  n  = 


sin(^)  _  cos(^) 


+l  +] 


The  skin  fnction  coefficients,  normalized  by  S  =  7d}^  and  c  =  2b,  are  yielded  by 
discretely  integrating  the  following  equations  in  a  cylindrical  coordinate  system: 


qa=^  1=1  ;=1 


=  2£2(^„C^rv)  JA^^,Ax,  /  2a 
1=1  >=1 


r  =Mj^,  Msp 
q^Sc’ 


2ZZ[(^=oC^at)  . (-  X,  /  2a •  a  +  z,  / 2a- 7)]a^,Ax,,  /  2a  B.14 
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Where  the  unit  surface  velocity  vector  is  given  by: 


cosi 


V  = 


w"  +  1 


/wsin(^)cos(;') 

V/w  +1 


+  cosi 


'Wsin(r) 


y  + 


wcosi 


(^)cos(r) 


+  sm 


+1 


Wsin(r) 


B.15 


The  surface  and  local  slope  of  a  prolate  spheroid  comes  from  the  following  relations; 


2  2 
X  r 


Prolate  Spheroid:  —  +  —7  =  1 


a 


1 

slope  m  =  —  =  - 
dx 


bx 


V 

a 

B.16 


Note:  The  forces  are  summed  over  half  the  spheroid,  ^  =  0  — >  1 80° ,  and  doubled. 

The  y-direction  forces  and  the  roll  and  yaw  moments  are  neglected  zero  due  to 
symmetry. 
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APPENDIX  C. 

MATLAB  PROGRAMS  TO  INTEGRATE  AERODYNAMIC  FORCES  OVER 
THE  SURFACE  OF  A  PROLATE  SPHEROID 


4m 


Page  1. 


%  This  Matlab  M-file  script  performs  a  first  order  (linear  approximation) 

%  integration  of  pressure  forces  over  a  6:1  prolate  spheroid.  Central  Differencing 

-  .  .  .  n  ^ _  I.V.  i  i  Data  iS  inOUt  frOITl 

% 

% 

% 


of  location  is  used  for  the  first  order  integration  routine.  Data  is  input  from 
AGARD  AR-303  Test  C-2  as  rotation  angle,  x/c,  and  Cp,  Data  is  for  two  test 
conditions,  AOA  =  10  and  29.7  degrees. 


Re  s  7.7xl0e6,  Vinf  =  55  m/s: 


clear 

%  Load  in  experimental  data  for  AOA  =10  degrees 
%  M  =  0.162 

fidl=fopen{ 'icpl Draws 'r+M ;  %open  file  for  prining  step  data  for  error  checking 

fprintf(fidl,'i  j  S  ^  S,  “ 

nx  ny  nz  ds  dN  dA  dm  \nM ; 

data  =  [cpl0data(:,2)  cplOdata ( : ,4)  cplOdata ( : , 5) ] ;  %  Extract  coluims  2,4,5. 
nphi  =  40;  nxc  =  42; 


m  =  0;  N  =  0;A  =  0;  Si=0; 
a  =  0.5;  b  =  0.5/6; 

S  =  pi*b^2; 
for  i  =  l:nphi 

for  j  =  l:nxc 

1  =  (i-l}*nxc+j; 
phi  =  data (1, 1) ; 

X  =  data{l,2)-a; 

Cp  =  data{l,3); 
r  =  b*sqrt(l-x'"2/a'"2)  ; 

M  =  -b*x/(a^2*sqrt(l-x^2/a""2+.000001))  ; 

2  =  ~r*cos(phi); 

nx  =  -M/sqrt(M^2+l) ; 

ny  =  sin (phi/57. 296) /sqrt{M^2+l) ; 

nz  =  -cos{phi/57.296)/sqrt(M'^2  +  l); 

nt=sqrt (nx^2+ny^2+nz^2) ; 

if  j  ==  1 

dx  =  data  (1+1, 2) /2 ; 
elseif  j  ==  nxc 

dx  =  (2*a-data{l-l,2))/2; 

else 

dx  = 

diate  pressure  ports 
end 


%Initialize  #  of  rotation  steps  and  pressure  ports 
%  Initialize  summed  forces  to  zero 
%  a  and  b  for  6:1  Prolate  Spheroid 
%  Reference  area  -  max  cross  section 


%  dx  at  first  pressure  port 
%  dx  at  last  pressure  port 


(data(l+l,2)-data(l-l,2) )/2;  %  central  differencing  at  interme 


if  i  ==  1 
dphi 
elseif  i 


data( (i*nxc+l) ,1) /2; 
==  nphi 


dphi  =  (180“data( ( (i-1) *nxc) ,1) )/2; 

dphi  =  (data ((i*nxc+l),l) -data (((i-l)*nxc),l}}/2; 


dA  dm]  ; 


else 
end 

dS  =  r*dphi/57.296*dx*sqrt(M^2+l) ; 
dN  =  2* (-Cp) *dS*nz; 
dA  =  2*  (-Cp) *dS*nx; 
dm  =  2* (-Cp) *dS* (-x*nz+z*nx) ; 

N  =  N+dN; 

A  =  A+dA; 
m  =  m+dm; 

raw(l,;)=(i  j  1  phi  dphi  x 
Si=Si+2*dS; 


Mr  dx  Cp  nx  ny  nz  dS  dN 


end 

end 

CN_J^OA10  =  N/S 
CA^OAIO  =  A/S 
CM^OAIO  =  m/(S*2*b) 

fprintf (fidl, '%3.0f  %3.0f  %5.0f  %5.2f  %6.4f  %4.3f  %5.2f 

%4.3f  %4.3f  %4.3f  %8.7f  %8,7f  %8.7f  %8.7f\n  ^,rawO; 

fprintf (fidl,' i  j  1  Pbi  dphi  ,  x  M 

nx  ny  nz  ds  dN  dA  dm  \n' )  ; 

fclose ( 'all ' ) ; 


%6.3f  %6.4f  %4.3f 


dx 


Cp 


! 
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<Jun  4109701:^^ _ tcf„pro{ate.in _ Page1 


%  This  Matlab  M-file  script  performs  a  first  order  (linear  approximation) 

%  integration  of  skin  friction  over  a  6:1  prolate  spheroid-  Central  Differencing 
%  of  location  is  used  for  the  first  order  integration  routine.  Data  is  input  from 
%  AGARD  AR-303  Test  C-2  as  rotation  angle^  x/c,  Cf  and  gamma  (crossflow  angle)  .  Dat 
a  is  for  two  test 
%  conditions,  AOA  =  10 
clear 

%  Load  in  experimental  da:ta  for  AOA  =  10  degrees.  Re  =  7,7xl0e6,  Vinf  =  55  m/s: 

%  M  =  0.162 
load  cflOreorder 
%fidl=fopen('icfl0raw',  'r+M  ; 

%fprintf  (fidl, 'i  j  1  phi  dphi  x  r  m  dx  Cf 

gamma  vx  vy  vz  dS  dN  dA  dm  \n')  ; 


data  =  [cflOreorder (: ,2) 
act  columns  2,4,5 
nphi  =  74;  nxc  =  12; 
Nsf=0;  Asf=0;  msf=0; 
m  =  0;  N  =  0;A  =  0;  Si=0; 
a  =  0.5;  b  -  0.5/6; 

S  =  pi*b'^2; 


cf lOreorder { : , 1)  cflOreorder (:, 3)  cflOreorder (;, 4) ] ;  %  Extr 

%Initialize  #  of  rotation  steps  and  pressure  ports 

%  Initialize  summed  forces  to  zero 
%  a  and  b  for  6:1  Prolate  Spheroid 
%  Reference  area  -  max  cross  section 


for  i  =  l:nxc 

for  j  =  l;nphi 

1  =  (i-1) *nphi+j ; 
phi  =  data (1,1); 

X  =  data (1,2) -a; 

Cf  =  data (1, 3) ; 
gamma  -  data(l,4); 
r  =  b*sqrt(l-x^2/a^2)  ; 

M  =  -b*x/ (a^2*sqrt(l“X^2/a^2+. 000001) ) ; 
z  =  -r*cos(phi); 

vx  =  cos (gamma/ 57.3) /sqrt (M^2+l)  ; 

vy  ss  M*sin {phi/57. 296) *cos (gamma/57. 3) /sqrt {M^2+l) +cos (phi/57.3) *sin(gamma/ 


vz  =  -M*cos (phi/57.296) *cos (gamma/57.3) /sqrt {M^2+l) +sin {phi/57. 3) *sin(gamma 


/57.3); 

nt=sqrt  (vx^2+vy'^2+vz'"2)  ; 
if  i  ==  1 

dx  =  data(nphi+l,2) /2;  %  dx  at  first  hot  film  sensor 

elseif  i  ==  nxc 

dx  =  (2*a-data( ( (i-1) *nphi) ,2) ) /2  ;  %  dx  at  last  hot  film 


sensor 

else 

dx  =  (data{ (i*nphi+l) ,2) -data ( { (i-1) *nphi) ,2) ) /2  ; 

end 

if  j  ==  1 

dphi  =  data ( (1+1), l)/2; 
elseif  j  ==  nphi 

dphi  =  (180-data(l-l,l))/2; 

else 

dphi  =  (datad  +  l,  1) -data (1-1,1)  )/2; 

end 

dS  =  r*dphi/57.296*dx*sqrt(M^2+l)  ; 
dN  =  2* (Cf) *dS*vz; 
dA  =  2* (Cf) *dS*vx; 
dm  =  2* (Cf ) *dS* (-x*vz+z*vx) ; 

Nsf  =  Nsf+dN; 

Asf  -  Asf+dA; 
msf  =  msf+dm; 

%  raw(l,:)  =  [i  j  1  phi  dphi  x  Mr  dx  Cf  gamma  vx  vy  vz 

dS  dN  dA  dm]  ; 
end 


end 

CNsf^OAlO  =  Nsf/S 
CAsf^OAlO  =  Asf/S 
CMsf_AOA10  =  msf/(S*2*b) 

%fprintf (fidl, '%3.0f  %3.0f  %5.0f  %5.2f  %6.4f  %4.3f  %5.2f  %6.3f  %6.4f  %6.5f 

%4.1f  %4-3f  %4.3f  %4.3f  %8.7f  %8.7f  %8.7f  %8.7f\n  SrawM; 

%fclose( 'all' ) ; 
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APPENDIX  D.  REPRESENTATIVE  CMARC/PMARC  SPEED  TEST  FILE 
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APPENDIX  E.  MATLAB  PROGRAM  FOR  REORDERING  AGARD  DATA  FILE 


97 
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Printed  by  pollard  trom  hawkeye 


APPENDIX  F.  CMARC  PROLATE  SPHEROID  INPUT  FILE 
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Printed  by  pollard  from  hawkoye 
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Printed  by  pollard  from  hawkeye 
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Printed  by  pollard  from  hawkeyo 
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Printed  by  pollard  from  hawkeye 
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Printed  by  pollard  from  hawkeye 
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Primed  by  pollard  from  hawkeye 


APPENDIX  G.  CMARC/PMARC  DATA  EXTRACTION  PROGRAM 


107 
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APPENDIX  H.  LOFTSMAN  INPUT  FILES 


^«n  1.3  m7.0&3& 


fogfu&a,ift 


Page.1 


BOX  MOLDLINES  DATA  TEMPLATE 

File  name:  fogfusa 
Last  revision:  4/12/97 


BOTTOM  WATERLINE 
Segments ;  3 


Fore  end 
Aft  end 
Comer 
Curvature 

0,6.5 

12,0  53.0,0  53.5,5.75 

0,0  S  53,5,0 

.69  0.95 

WAIST  WATERLINE 

Segments : 

1 

Fore  end 

0,6.5 

Aft  end 

53.5,6.5 

Comer 

S 

Curvature 

TOP  WATERLINE 

Segments : 

7 

Fore  end 

0,6. 

5 

Aft  end 

8,9. 

3  15.2,9.3  21.6,13.0  29.0,14.6 

44.6,11.6 

53.0,11.5  53.5,6. 

Comer 

0,9. 

3  S  S  24.3,14.6 

35,4,14.6 

S  53.5,11. 

Curvature 

0.7 

0.71 

0.81 

.95 

MAXIMUM  BUTTLINE  DISTANCE  FROM  PLANE  OF  SYMMETRY 

Segments : 

6 

Fore  end 

0,0 

Aft  end 

1,3  1,3  22,4.5  43.6,4.5 

53.1,1,3 

53.5,0 

Comer 

0,2.9  S  S  S 

S 

53.5,1.3 

Curvature 

.9 

0.8 

0.95 

BOTTOM  K  FACTOR 

Segments :  3 

Fore  end 
Aft  end 
Comer 
Curvature 


0,0.93  • 
12.0,0.98 
S 


43.6,0.98 

S 


53.5,0.95 

S 


TOP  K  FACTOR 

Segments :  4 

Fore  end 

Aft  end 

Corner 

Curvature 

0,0,90 

15.20,0.95 

S 

24,1.0 

S 

44.65,1.0 

S 

53.5,0.95 

S 

BUTTLINE  AT 

PLANE  OF  SYMMETRY 

Segments :  0 


109 


Pag&1 


iuni:i  myoscg^  ■ _ fogwlnga.wi 

NPS  FROG  UAV  Main  Wing  -  Loftsman  Input  File 
Date:  S/29/91 
Breaks :  5 
Break  1 

Axis:  24.65,0,13.1 
Axis /chord:  0 
Chord:  20.0 
Incidence:  4.5 
Cant:  0 

Section  file:  N2415 
T/C  ratio:  0.1500 
Spars :  0 

Panel  rib  angles:  0,999.0000,0.0000 
Break  2 

Axis:  24.65,6,13.1 
Axis /chord:  0 
Chord:  20.0 
Incidence:  4.5 
Cant :  0 

Section  file;  N2415 
T/C  ratio:  0.1500 
Spars :  0 

Panel  rib  angles;  0,999.0000,0.0000 
Break  3 

Axis:  24.65,31.5,13.1 
Axis /chord;  0 
Chord;  20.0 
Incidence:  4.5 
Cant :  0 

Section  file:  N2415 
T/C  ratio:  0.1500 
Spars :  0 

Panel  rib  angles:  0,999.0000,0.0000 
Break  4 

Axis:  24.65,53.0,13.1 
Axis /chord:  0 
Chord:  20.0 
Incidence;  4.5 
Cant :  0 

Section  file:  N2415 
T/C  ratio:  0.1500 
Spars :  0 

Panel  rib  angles:  0,999.0000,0.0000 
Break  5 

Axis:  24.65,61.0,13.1 
Axis/chord:  0 
Chord:  18.5 
Incidence:  4.5 
Cant :  0 

Section  file:  N2415 
T/C  ratio;  0.1500 
Spars ;  0 


Junti  1S&7  05:24 _ lQg|eftpOd.lft _ Pagel 

FROG  UAV  ENGINE  NACELLE 


File  name;  fogenpod 
Last  revision:  4/13/97 


BOTTOM  WATERLINE 

Segments :  4 

Fore  end  16.5,20.4 

Aft  end  18.2,18.6 

Comer  16.6,19.6 

Curvature  0.79 

21.0,16.8 

19.15,17.35 

0.83 

31.0,15.75 

23,8,15.9 

0.72 

43.0,  16.8 

35,6,15.65 

0.73 

WAIST  WATERLINE 

Segments :  1 

Fore  end 

16.5,20.4 

Aft  end 

43.0,16.8 

Comer 

S 

Curvature 

TOP  WATERLINE 

Segments ;  4 

Fore  end 

16.3,20.4 

Aft  end 

18.45,22.1 

27.0,21.75 

35.0,19,8 

43.0,16.8 

Comer 

16.75,21.3 

21.4,22.5 

30.4,21.3 

38.3,18.75 

Curvature 

0.79 

0.80 

0.70 

0,75 

MAXIMUM  BUTTLINE  DISTANCE  FROM  PLANE  OF 

SYMMETRY 

Segments :  4 

Fore  end 

16.5,0 

Aft  end 

18.2,1.6 

23.0,2.3 

40.8,2.3 

43.0,0 

Comer 

16.5,0.70 

20.1,2.25 

S 

43.0,2.3 

Curvature 

0.72 

0.75 

0.90 

BOTTOM  K  FACTOR 

Segments :  4 

Fore  end 

16.5,0.707 

Aft  end 

18.2,0.707 

24.0,0.93 

42.0,0.93 

43.0,0.75 

Comer 

S 

20,0.93 

S 

S 

Curvature 

0.9 

TOP  K  FACTOR 

Segments ;  4 

Fore  end 

16.5,0.707 

Aft  end 

18.45,0.707 

24.5,0.93  42.0,0.93 

43.0,0.75 

Comer 

S 

20.3,0, 

93  S 

S 

Curvature 

0.9 

BUTTLINE  AT  PLANE  OF  SYMMETRY 
Segments :  0 
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jun  11 1^7  05:27  Jft _ Pagei 


FROG  UAV  ENGINE  PYLON  (Lofted  as  A- Body  Type) 

“Basic  pylon  model  modified  so  as  not  to  have  a  top  and  bottom 
“Single  strip  which  is  the  side  of  the  pylon. 

File  name:  FOGPYLOl 
Last  revision:  4/23/97 


Strips ;  1 
Sym:  Y 

MIB 


Segments:  1 

Fore  end 

25.8,0 

Aft  end 

37.7,0 

Comer 

25.8,3.8 

K  factor 

0.71 

MIW 

Segments :  1 

Fore  end 

25.8,14.31 

Aft  end 

37.7,13.65 

Comer 

31.0,15.45 

K  factor 

0.72 

CIB 

Segments :  S 

CIW 

Segments:  S 

K1 

Segments ;  S 

M2B 

Segments :  ^MIB 

M2W 

Segments :  1 

Fore  end 

25.8,16.01 

Aft  end 

37.7,16.1 

Corner 

33.65,15.35 

K  factor 

0.65 
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JUntI  1Q&7  05:22 _ fogbOQmJft 

FROG  UAV  Tail  Boom 


File  name:  frogboom 
Last  revision:  4/28/97 

4/28:  added  rounded  start  and  finish  to  close  ends 

BOTTOM  WATERLINE 


Segments :  3 

Fore  end 
Aft  end 
Comer 
Cuirvature 


53.5,9.375 

54.8.5 

53.5.8.5 
0.707 


WAIST  WATERLINE 


88,8.5 

S 


Segments :  1 

Fore  end 
Aft  end 
Comer 
Curvature 


53.5,9.375 
88.5,  9.375 
S 


TOP  WATERLINE 
Segments :  3 


Fore  end 
Aft  end 
Comer 
Curvature 


53.5,9.375 

54.0,10,25 

53.5,10.25 

0.707 


88,10.25 

S 


MAXIMUM  BUTTLINE  DISTANCE  FROM  PLANE  OF  SYMMETRY 


Segments :  3 

Fore  end  53.5,0 

Aft  end  54,0.875 

Comer  53.5,0,875 

Curvature  0.707 

BOTTOM  K  FACTOR 

Segments:  1 

Fore  end  53.5,0.707 

Aft  end  88.5,  0.707 

Comer  S 

Curvature 

TOP  K  FACTOR 

Segments:  1 

Fore  end  53,5,0.707 

Aft  end  88.5,  0.707 

Corner  S 

Curvature 

BUTTLINE  AT  PLANE  OF  SYMMETRY 
Segments;  0 


88,0.875 

S 


Page1 


88,5,9.375 

88.5,8.5 

0.707 


88.5,9,375 

88,5,10.25 

0.707 


88,5,0 

88,5,0.875 

0.707 
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05:26 _ ,  fOghtaH.wl  _ Pgg&l 

FROG  Horizontal  Tail 
Date:  4/14/97 
Breaks :  2 
Break  1 

Axis:  82.5,0,8.09 
Axis/chord:  0 
Chord:  13.5 
Incidence:  0 
Cant ;  0 

Section  file;  N0006 
T/C  ratio;  0.06 
Spars :  0 

Panel  rib  angles;  0,999.0000,0.0000 
Break  2 

Axis:  86.5,19.875,8.09 
Axis /chord:  0 
Chord:  9.55 
Incidence:  0 
Cant :  0 

Section  file:  N0006 
T/C  ratio :  0.06 
Spars :  0 

Panel  rib  angles;  0,999.0000,0.0000 
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<^un  n  iwe&ge  fo^vert>wi _ 


FROG  UAV  Vertical  Tail  -  LOFTSMAN  input  file 
Date:  AflA/91 
Breaks :  2 
Break  1 

Axis:  77.5,0,10.4 
Axis /chord:  0 
Chord:  20 
Incidence:  0 
Cant:  90 

Section  file:  N0006 
T/C  ratio:  0.06 
Spars :  0 

Panel  rib  angles:  90,0,999 
Break  2 

Axis:  92.35,0,25.15 
Axis /chord:  0 
Chord:  10 
Incidence:  0 
Cant:  90 

Section  file:  N0006 
T/C  ratio;  0.06 
Spars :  0 

Panel  rib  angles:  90,0,999 
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Printed  by  pollard  from  hawkeye 


APPENDIX  L  FROG  UAV  CMARC  INPUT  FILE 
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Printed  by  pollard  from  hawkeye 


119 


Printed  by  pollard  from  hawkeye 
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121 


122 


Printed  by  pollard  from  hawkeye 
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11.4914 
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Printed  by  pollard  trom  hawkeye 


126 


Printed  by  pollard  from  hawkey© 


A 


Printed  by  pollard  from  hawkeye 


L 


Printed  by  pollard  from  hawkeye 
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i 


Printed  by  pollard  from  hawkeye 
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Printed  by  pollard  from  hawkeye 
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i 


Printed  by  pollard  from  hawkeye 
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